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ABSTRACT 


Automated Guided Vehicles (AGVs) use different techniques to help locate their 
position with respect to a point of origin. This thesis compares two approaches that 
utilize a binary track laid on the floor for the AGV to follow. Both approaches use 
equally spaced n-tuples on the track that the AGV can use to compute its position. Both 
approaches also have the special feature that every n-tuple on the binary track is unique 
and can be used to designate the position of an AGV. The first approach, developed by 
E.M. Petriu, uses a Pseudo-Random Binary Sequence (PRBS) as a model for the binary 
track. In the second approach, we use a Greedy DeBruijn Sequence (GDBS) as a model 
for the binary track. Unlike the PRBS model, the GDBS model has a natural ordering 
which can be used to determine the position of the AGV more quickly and efficiently 


than the PRBS model. 
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EXECUTIVE SUMMARY 


The greedy DeBruijn sequence (GDBS) of span n is a binary string of length 
2” where every n-tuple in the sequence is unique. Originally the sequence was defined 
by the rule: Start with n-1 zeros and then add a 1. Continue adding 1’s as long as the n- 
tuple created has not been seen before in the sequence. If it has already appeared, add a 
zero instead. Continue by adding Is. If neither a 0 nor a 1 can be added, then stop. The 
process will generate a sequence of length 2” [Mar]. The sequence can also be 
constructed more efficiently by concatenating a list of combinatorial objects called 
necklaces in lexicographic order. [Fr] The motivation behind this project is to try to 
determine a mapping between the set of necklaces for the greedy DeBruijn sequence and 
the natural code, much like the mapping that exists between the Gray code strings and the 
natural code. Finding this necklace-to-natural code mapping would enable us to 
determine the absolute position of any arbitrary n-tuple in a greedy DeBruijn sequence. 
Although we do not solve the mathematical problem in determining this mapping, we are 
able to apply the properties of this DeBruijn sequence to construct a very efficient 
implementation in helping an Automated Guided Vehicle (AGV) determine its position 
on a long DeBrujjn track. Since every n-tuple in the DeBruijn sequence is unique, the n- 
tuples can be used to designate unique positions for the AGV. 

In this thesis, we compare two approaches to the AGV implementation. E.M. 
Petriu developed the first approach. He uses a Pseudo-Random Binary Sequence (PRBS) 
generated by a linear shift register as the basis of his binary track. He implements a 
series of n-tuples, called milestones, on the track that provide the AGV with information 
about its location. We advocate the use of the GDBS as the binary track since it has 
certain advantages over the PRBS. Like Petriu, we use certain n-tuples, called signposts, 
to give the AGV information about its position. However, the way the signposts are 
generated and they way they are used differs significantly from Petriu’s approach. With 
the GDBS, the AGV can determine its location in a far more efficient manner than can 


the method used by Petriu. 


XV 


Chapter II gives the background on the mathematical structures of necklaces. 
There the Necklace Algorithm developed by Fredricksen, to generate necklaces 
efficiently, is introduced. [Fr] This algorithm plays an important role in helping the AGV 
determine its absolute position. There the GDBS is constructed from the necklaces as 
generated by the Necklace Algorithm. 

Chapter III introduces Petriu’s method. It describes the sequential-parallel 
approach Petriu uses to determine the AGV’s absolute position. We also discuss the 
equipment and time costs associated with his scheme and the computation needed to 
determine the optimum number of milestones that will be used. 

Chapter IV describes our method and the results obtained. We introduce the 
initial steps the AGV must take to determine its location. We describe how signposts are 
generated and where they are placed on the GDBS. We compare the number of 
milestones versus the number of signposts and find they are similar in both methods but 
the GDBS method has computational advantages over Petriu’s PRBS scheme. The 
GDBS scheme can also be used to potentially optimize the distribution of more than one 
AGV on the GDBS track. Finally, we describe an alternate means of signpost generation 
and compare it to the approach we decided to use. 

Chapter V describes current problems of a mathematical nature that 
simultaneously give insight into necklaces and their structure and significantly improve 
the AGV implementation. By solving the necklace-to-natural code-mapping problem, the 
AGV would not even need an external binary track to follow. There are two fruitful 
approaches that were initiated to try to solve this problem but more work needs to be 
done. 

Appendix A provides a background on linear shift registers and how they generate 
the PBRS. Appendix B contains the code used to execute the Necklace Algorithm, 
generate necklaces and signposts, and gather statistics needed to analyze the structure of 


necklaces. 
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I. WHAT IS A NECKLACE? 


A. BACKGROUND 

This chapter provides background material on using a greedy DeBruijn sequence, 
instead of a PBRS, as a bit track for the Automated Guided Vehicle (AGV). This 
DeBruijn sequence, generated via the Necklace Algorithm, has some properties that give 
it an advantage over the PBRS in helping determine the position of an arbitrary n-tuple of 


that sequence. 
B. STRUCTURE OF NECKLACES 


Let A, represent the set of all n-bit binary strings. The total number of possible 
strings in A,is 2". We define a mapping ¢:A, — B, where the set B, is the set of all 
necklaces over A,. By anecklace, we mean a linear representation of the collection of 


all circular permutations of a binary n-tuple into another binary n-tuple. For example, for 


n=5, let a,= 01101 be an arbitrary binary string in A,. Then the following binary 


strings are all left circular permutations of a,: 


Figure 1. | Circular Permutations of 01101 


Note that each binary string is a result of cyclically rotating the previous string 
one position to the left. Also, note that the last cyclic permutation just reproduces the 
original string a,. From this, we can see that each binary string in this list is some 


number of cyclical shifts of any other binary string in the list. Although each binary 
1 


string is distinct in A,, we can represent this whole list of binary 5-tuples as an 
equivalence class (the elements of B, ) of some binary string in the list. We pick a 


particular binary string to be the linear representation of this equivalence class. We call 


this representative a necklace, and we define it to be the element of A, that has the 


largest decimal value. (See Table 1.) 

















Binary String | Decimal Equivalent 
01101, 13 
11010. 26 
10101. 21 
01011, 11 
101102 22 














Table 1. | Decimal Equivalents of Elements of Equivalence Class for 01101 


We see that 110102 has the greatest decimal equivalent so this is the “necklace” 
for this equivalence class. Note that the necklace always ends in a “0” (except for the 
necklace of all “1”’s) since a “1” at the end can always be cyclically shifted to the front of 
the binary string to give a larger decimal equivalent. Counting from the Most Significant 
Bit (MSB), or the leftmost bit, a necklace representative contains its greatest number of 
leading ones before the first zero. 


The crucial point here is that every possible binary string in A, (for every n) will 


be a member of some necklace’s equivalence class. Further, it is possible to produce a 
list of necklaces whose equivalence classes are disjoint (no two distinct equivalence 
classes can contain any binary string in common). One obvious way to tell if two 
necklaces are from two different equivalence classes is to see if they have different 
densities (different number of ones). The certain way to ensure that two equivalence 
classes are disjoint is to test if one necklace can be circularly permuted into the other one. 


For example, let b,= 110010 and b,= 110100 be two necklaces for n = 6. Neither can be 


obtained from the other by circular permutation although both have the same density. 


Both have different equivalence classes as seen in Table 2: 




















Equivalence Class for 110010 | Equivalence Class for 110100 
110010 110100 
100101 101001 
001011 010011 
010110 100110 
101100 001101 
011001 011010 














Table 2. | Comparison of Equivalence Classes 


One could randomly choose binary strings in A, and search for necklaces, but this 


would be an inefficient procedure at best (especially for large n). Fortunately, an elegant 


algorithm has been developed to generate all the necklaces for any n [Fr]. 


c. NECKLACE ALGORITHM 


The Necklace Algorithm [Fr] generates the list of all necklaces in lexicographic order 


from the all ones necklace to the all zeros necklace. We use a subroutine, called the 


©—Algorithm, to generate a binary string, which may be a necklace. 


©-Algorithm [Fr]: 


1. 


Given an arbitrary binary sequence of length n (where the indices ascend 
from left to right), read the binary string from right to left and find the 
largest index corresponding to the bit that is equal to “1” and any 
remaining subsequent bits (bits of a higher index) are all Os. Denote this 
largest index by j. Ifj is equal to n, there will be no zeros following it. 


Denote this binary n-tuple as the parent string. Copy the contents of this 
parent string into another string called the child string. Subtract “1” from 
the bit corresponding to the j” index in the child string. This bit will 
become a zero. Discard the remaining n-j bits (all zeros) of higher index 
following this bit. 


Since we discarded the remaining n-j bits, we may need to fill them with 
bits so that the child string has a length of n. We do this by copying n-j 


bits from the beginning of the child string into the remaining n-j “spaces” 
as many times as necessary to complete an n-long string. 


The following theorem [Fr] states: 
If a=4a,4,,...,a, is a necklace, # 000...0, then one of O(a), 0’ (a),....0" (a) is a 


necklace. 

As an example, let the 11110 be the parent string. We copy the contents of this parent 
string into a child string. Then, 7 = 4. Subtracting “1” from the bit corresponding to this 
index, we have 

11110 (child string) 
za | 
1110_ (child string). 
To fill the remaining position, we need to copy one bit from the beginning of the child 


string so we have: 


11101. 
Necklace Algorithm [Fr]: 
1. Start with the n-bit binary string of all ones denoted by 1” . This is the first 


necklace and the first binary parent string. 


2: To find the (+1)” necklace fori >1, apply the ©—algorithm to the 
i” necklace and denote the execution of this algorithm as 6’. 

3: If the j” index of the parent string divides n, the child string is a necklace. 
Then increment i and return to step 2. Otherwise it is nota necklace. The 


child string now becomes the parent string for the next execution of the O- 
Algorithm. 


4. If j does not divide n, then apply 0’, @’,...,0°, where k denotes the least 
number of times needed to apply the @—algorithm to the i” necklace until j 
divides n. The resulting child string is a necklace so increment i. 


a. If the all zeros necklace 0” has not appeared in steps 2,3 and 4 above, 
return to step 2. 


The number of bad n-tuples (or non-necklaces), B, , between any pair of consecutive 


necklaces generated by the Necklace Algorithm is bounded by B, < EG - D Z,, where 


Z, == 52 g(a) (1.1) 


d\n 


Z,,1s the number of necklaces, d is a divisor of n, and g(d) is Euler’s totient function 


which represents the number of integers that are relatively prime to d, including “1” [Go]. 
By the Necklace Algorithm, we generate a lexicographic list of necklaces whose 
equivalence classes together contain all possible 2” binary strings of length n. The 
following example demonstrates the Necklace Algorithm. 
Let n = 6. Starting with the initial necklace and parent string: 
111111 
note that the index j = 6 and there are no zeros following the last bit. 
Also, the index of this parent string (j = 6) divides 6, so the child formed from the all 
ones parent string will also be a necklace. Subtracting 1 from the last bit we have: 
111110 
Since the string has 6 bits, the ©—Algorithm need not fill in any missing bits. 
This is the second necklace and the new parent string. Note that the index of the last non- 
zero bit is 5. This index does not divide 6 so the child string formed from the parent 
string 111110 is not a necklace. The ©—Algorithm procedure yields: 
111110 (parent string) 
_-1 
11110_ (child string) 


To fill the remaining position, we copy one bit from the beginning of the child string so 


| 4 


111101 


we have: 


To see how this child string is not a necklace, just cyclically shift the last bit to the 


beginning of the string to produce 


111110 
Note that this is the parent string that was a necklace. Remember that a necklace is the 
representative of the equivalence class that has the largest decimal value. It is easily seen 
that 111101 belongs to the equivalence class of 111110 and that 111110.>1111012. To 
continue, note that for the string 111101, 7 = 6. This obviously divides six so the child 
string of 111101 is a necklace although the string 111101 is not a necklace. To see the 
necklace, subtract “1” from 111101 to produce 
111100 
which is clearly a necklace. (Further this necklace cannot be cyclically rotated to produce 
111110. So the equivalence class of 111100 is distinct from the equivalence class of 
111110). 
Continuing in the same manner, we generate a list containing both binary 
necklaces and non-necklaces from 111111 to 00000. 
111111 - necklace (j” index = 6) 
111110 - necklace (j” index = 5) 
111101 - non-necklace (j" index = 6) 
111100 - necklace i" index = 4) 
111011 - non-necklace (j” index = 6) 
111010 - necklace (j” index = 5) 
111001 - non-necklace (j” index = 6) 
111000 - necklace (j” index = 3) 
110110 - necklace (j index = 5) 
110101 — non-necklace gi" index = 6) 
110100 - necklace i" index = 4) 
110011 - non-necklace (j" index = 6) 
110010 - necklace (j” index = 5) 
110001 - non-necklace (j” index = 6) 
110000 - necklace (j” index = 2) 
101010 — necklace (j” index = 5) 
101001 - non-necklace (j" index = 6) 


101000 - necklace (j index = 3) 

100100 - necklace (j” index = 4) 

100010 - non-necklace (j" index = 5) 

100001 - non-necklace i" index = 6) 

100000 - necklace (j” index = 1) 

000000 - necklace (end of algorithm) 

With the necklace algorithm, we only need to generate and search through 23 (instead of 
2° = 64 ) binary strings to find out which ones were necklaces. Of these 23, exactly 14 
are necklaces. 

D. DEBRUIJN SEQUENCES 


A binary DeBruijn sequence of span n contains every possible n-tuple of length n 


in a cycle of length 2”. The necklace x,x,...x, (where eachx, € {0.1}) has a period of 
Tength:.2 18 4X) X= HH wk, Hips hyn OO 
NyXq ee Xq = gag Bgaa-Xog = Xo qsj%oG49+-%gq =- 18 the periodic sub-necklace of periodd . 


Lyndon words are necklaces that exhibit no sub-period behaviour or necklaces with 
multiple cycles removed. We are only interested in concatenating all the Lyndon words 
from the Necklace Algorithm, in lexicographic order, to generate a DeBruijn sequence. 
The length of a Lyndon word is equal to one period of a necklace. Table 3 lists all the 
necklaces along with their periods and respective Lyndon words for the case where n = 6. 
(We list those Lyndon words in bold face that are shortened by virtue of being periodic 


necklaces). 























Necklace | Lyndon Word | Period 
111111 1 1 
111110 111110 6 
111100 111100 6 
111010 111010 6 
111000 111000 6 
110110 110 3 






































Necklace | Lyndon Word | Period 
110100 110100 6 
110010 110010 6 
110000 110000 6 
101010 10 2 
101000 101000 6 
100100 100 3 
100000 100000 6 
000000 0 1 

















Table 3. | Necklaces and Lyndon Words for n = 6 


To determine the period of a necklace, simply count the number of bits until the 
pattern of ones and zeros begins to repeat. It is easy to see this process if one 
concatenates copies of the necklace with itself (where the carat “” operator is used for 
concatenation) and count how many bits are needed before the same pattern of zeros and 


ones is repeated. In the following examples 


L111114111111 =111111111111 
000000000000 = 000000000000 


the pattern repeats itself after only one bit so the period is one. The following sequence 


1101104110110 = 110110110110 


repeats the patterns of ones and zeros after three bits. So the period is three. However, 


the following sequence 


111100111100 = 111100111100 


requires six bits before the pattern repeats itself. So its period is six. We can now 
concatenate in order the Lyndon words found by the Necklace Algorithm to find the 


DeBruijn sequence. For n = 6 the DeBruijn sequence is: 


14111110%111100%111010%111000%110%110100%110010%110000%10%101000%100*%10000 
0“0 


The resulting sequence is 2° = 64 bits in length. Reading from left to right, the first 
string is 111111. Cycling one bit to the right and reading six bits yields the next string 
111110. Continuing in this manner until the last bit is reached yields all possible binary 
strings of length six that are contained in the above sequence exactly once. The binary 
strings 000001, 000011, 000111, 001111, and 011111 are read from the sequence as we 
cycle one bit at a time to the right from 000000 to 111111. 

The reason that we concatenate the Lyndon words, and not the necklaces 
themselves, is to ensure that every possible binary sequence of length six occurs only 
once. To see this, imagine if we concatenated the entire necklaces and not just the 


Lyndon words: 


1111114111110%111100%111010%111000%110110%110100%110010%110000 * 101010 
* 101000*100100%100000*000000 = 


1111111111101111001110101110001101101101001100101100001010101010001001001 
00000000000 (84 bits) 


If we read the sequence formed by concatenation from left to right, we see that the binary 
strings 111111 and 000000 occur multiple times. The binary strings 110110, 101010 and 
100100 occur multiply. By using the Lyndon words, every possible binary string of length 
six occurs exactly once in the DeBruijn sequence for n = 6. This is the same sequence as 


that formed by the greedy algorithm described in Chapter I. 


Note that n = 6 is a composite number (it has more factors other than one and itself) 
and that the factors of six (one, two, three and six) correspond to the periods of the 
necklaces (or the length of the Lyndon words). This is always the case for any n20. Ifn 
is prime (itself and 1 are the only factors), then the only necklaces that have periods less 
than n are the all 1s and the all Os necklace. All other necklaces have periods equal to n. 


Table 4 lists the necklaces for n = 5: 























11111 1 1 
11110 11110 i) 
11100 11100 5 
11010 11010 3) 
11000 11000 =) 
10100 10100 iS) 
10000 10000 =) 
00000 0 1 

















Table 4. | Necklaces and Lyndon Words for n = 5 


Concatenating only the Lyndon words 
1411110%11100%11010%11000%10100%10000%0 = 
11111011100110101100010100100000 
is the DeBruijn sequence of length 2° = 32 bits. 
E. ENUMERATION OF NECKLACES 

We organize the necklaces into classes according to their class number. Reading 
the necklace from left to right, the class number corresponds to the number of ones 
preceeding the first zero. When we list the necklaces for n = 6, they are listed in a 
decreasing order from 1111112 to 0000002. This is one of the special features produced 
by the necklace algorithm that we will use later on. Table 5 contains the necklaces, their 


class numbers and decimal equivalents for n = 6: 
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Necklace | Class Number | Decimal Equivalent 
111111 6 63 
111110 5 62 
111100 4 60 
111010 3 58 
111000 3 56 
110110 2 54 
110100 2 52 
110010 2 50 














101010 1 42 
101000 1 40 
100100 1 36 
100000 1 32 
000000 0 0 

















Table 5. | Necklaces, Class Numbers and Decimal Equivalents for n = 6 


For purposes of this presentation, we do not shorten the necklaces to their 
respective Lyndon words for those necklaces whose periods are less than n since we want 
the true decimal equivalence of the binary number. Note that the decimal equivalents 
(except for the first number) decrease by two until we reach the string 110000. The 
decimal equivalents of 46 and 44 are “missing”. The reason 46 and 44 are missing is that 
their binary equivalents are 1011102 and 1011002 respectively, and these strings are not 
necklaces since they can be cyclically rotated to produce 111010 and 110010 
respectively. As seen in the Table, the necklaces 111010 and 110010 have already 
appeared in the list earlier. So the strings 101110 and 101100 belong to the respective 
equivalence classes of 111010 and 110010. Recall that the list of necklaces will contain 
only one representative (i.e., the necklace) from each distinct equivalent class. All 
decimal equivalents that are “missing” as we continue on correspond to binary strings 


that are in equivalence classes of necklaces that have already appeared on the list. At 
11 


times, we may have more than one even decimal “missing” between two adjacent 
necklaces. We call this group of consecutive missing even decimals a “clump”. 

Define the threshold as the necklace that has the smallest decimal equivalent from 
the beginning of the list of necklaces and for which no “missing” has yet been 


encountered. This necklace is designated by p ones followed by all zeros where 
n—-l 
=| ——_ 1.2 
P 5 (1.2) 
6 


Forn=6, p= a =2 so our threshold is 110000 (highlighted in the Table above). 


The decimal equivalents have no missing values until the class number of the necklace is 
strictly less than p. Any necklace whose class number is strictly less than p (i.e., one 
and zero for n = 6) will have “missing” decimal equivalents. 

Define cardinality of a class as the number of necklaces in a given class. For n = 
6, the cardinality of class one is four. The entries in the Table below are the cardinalities 
for a given class number k = n-q and a given necklace length n and an index g. The totals 
represent the total number of necklaces generated by the necklace algorithm for a given 


n. For n=6, the class numbers and their cardinalities are given in Table 6. 





ClassNo. |6/5/4/3/2/1/0 
Cardinality | 1) 1|1)2|)4/4]1 
































Table 6. | Class Number vs. Cardinality for n = 6 


Summing the cardinalities produces 14 necklaces for n = 6, which, not surprisingly, 
corresponds to the number of necklaces generated by the Necklace Algorithm. We can 
also calculate the number of necklaces of length n without running the Necklace 
Algorithm by using equation 1.1 on page 5 to calculate Z,. Although the number of 


necklaces grows with increasing n, it is small compared to the number of total binary 


n 


p : pe hoes 
strings 2” since Z, ~ — [Go]. In other words, we can represent the entire binary space 
n 


1 


aM deb va. : ee 
containing 2" strings with —th of its members (represented as the % entries in Table 7 
n 


and Table 8). This difference between Z, and 2" is also displayed in Figure 2. 







































































q n=3 | n=5 | n=7 | n=9 | n=1l | n=13 | n=I5 | n=17 n=19 n=21 
0 1 1 1 1 1 1 1 1 1 1 

1 1 1 1 1 1 1 1 1 1 1 

2 1 1 1 1 1 1 1 1 1 1 

3 1 2 2 2 2 2 2 2 2 2 

4 2 4 4 4 4 4 4 4 4 

5 1 6 8 8 8 8 8 8 8 

6 4 14 16 16 16 16 16 16 

7 1 19 30 32 32 32 32 32 

8 9 50 62 64 64 64 64 

9 1 56 114 126 128 128 128 
10 18 178 242 254 256 256 
11 1 172 | 433 498 510 512 
12 40 635 944 1010 1022 
13 1 533 1640 1968 2034 
14 93 2262 3682 4016 
15 1 1646 6222 77177 
16 210 8072 14363 
17 1 5126 23610 
18 492 28828 
19 1 16035 
20 1169 
21 1 

Total | 4 8 20 | 60 | 188 | 632 | 2192 | 7712 | 27596 | 99880 

a 8 32 | 128 | 512 | 2048 | 8192 | 32768 | 131072 | 524288 | 2097152 
% 50 | 25 | 15.6 | 11.7 | 9.18 | 7.71 | 6.69 5.88 5.26 4.76 









































Table 7. Necklace Enumeration for Odd n 
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n=4 | n=6 | n=8 | n=10 | n=12 | n=14 | n=16 | n=18 n=20 n=22 
0 1 1 1 1 1 1 1 1 1 1 
1 1 1 1 1 1 1 1 1 1 1 
2 1 1 1 1 1 1 1 1 1 1 
3 2 2 2 2 2 2 2 2 2 2 
4 1 4 4 4 4 4 4 4 4 4 
=) 4 8 8 8 8 8 8 8 8 
6 1 11 16 16 16 16 16 16 16 
7 7 27 32 32 32 32 32 32 
8 1 33 59 64 64 64 64 64 
9 14 96 123 128 128 128 128 
10 1 101 223 251 256 256 256 
11 30 338 479 507 512 512 
12 1 305 844 991 1019 1024 
13 63 1202 1867 2015 2043 
14 1 940 3199 3914 4063 
15 142 4281 7276 8010 
16 1 2915 12128 15462 
17 328 15267 28374 
18 1 9078 46005 
19 765 54511 
20 1 28418 
21 1810 
22 1 

Total | 6 14 32 | 108 | 352 | 1182 | 4116 | 14602 52488 190746 
16 64 | 256 | 1024 | 4096 | 16384 | 65536 | 262144 | 1048576 | 4194304 
% 37.5 | 21.9 | 12.5 | 10.5 | 8.59 | 7.21 6.28 5.57 5.01 4.55 
Table 8. | Necklace Enumeration for Even n 
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Figure 2. Z, vs. 2” 
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YW. PETRIU’S METHOD 


A. BACKGROUND 

In this chapter we describe Petriu’s method for computing the Automated Guided 
Vehicle’s (AGV) position from a point of origin on a binary track. We also describe the 
hardware and software efficiencies of his method as well as its limitations. Petriu’s 
scheme only requires a 1-bit wide code track to be physically laid on the floor. This is an 
attractive alternative for absolute position measurement that requires a very long track. 
This method is used by high-resolution shaft encoders and Automated Guided Vehicles 
(AGVs) [Pe2]. He uses a pseudorandom binary sequence (PRBS) as the basis for the 
track that the AGV will follow on the floor. The PBRS is a maximal length linear 
sequence of length 2” —1 generated by a linear shift register (see Appendix A for more 
details). If we scan every possible n-tuple in this sequence with a window of length n, we 
see that every possible n-tuple (except the all 0 n-tuple) appears uniquely. If the AGV’s 
position is designated by an n-tuple, then it can use this uniqueness property to determine 
its exact position on this track [Pe2]. 
B. ABSOLUTE POSITION MEASUREMENT 

There are different methods to convert these pseudorandom n-tuples into a natural 
code (the binary equivalent of the exact distance from a point of origin). One conversion 
method uses a strictly parallel solution by implementing a ROM-stored lookup Table. 
This can increase hardware costs for large values of n [Pe2]. At the other extreme is a 
strictly serial solution by using a reverse feedback shift register to count the number of 
reverse feedback shifts it takes to get to the “zero position” code pattern. This is 
equivalent to counting the number of bits the n-tuple is from the origin. Although the 
hardware costs are not as high, it becomes prohibitively time consuming as n gets large 
[Pel]. Petriu uses a combination of the parallel and serial methods to achieve a more 
economic approach. This approach implements a set of evenly spaced n-tuples, called 
milestones, where each milestone and its distance from the point of origin (designated to 
be the MSB of the PBRS), is paired with a binary equivalent of the distance of the 
milestone’s MSB from the point of origin [Pe2]. 


1b 


In the PBRS, we have a total of 2” —1 positions, or n-tuples, that the AGV can 
occupy (the PBRS does not include the n-tuple consisting of all zeros). See Appendix A 
for further details. Since the milestones are uniformly distributed with a period of ¢ bits, 
the total number of milestones needed is 

w=|(2"*)/t | (2.1) 
Let p=m*t-+r designate the position of any n-tuple where m*t represents the position 
of the nearest “down the track” milestone, Q(m), and r represents the relative distance 
between this milestone and the n-tuple representing the AGV’s initial position [Pe2]. 

Petriu uses a sequential algorithm for the code conversion of the relative distance 
rand a parallel code conversion method for the milestone position m*t. The sequential 
algorithm counts the number of steps needed for the AGV to move from its initial 
position toward the origin (p=0) until it reaches the MSB of a milestone. The parallel 
code conversion scheme compares each unique n-tuple encountered during this “back 
stepping” against all possible milestones to see if the binary patterns match [Pe2]. The 
parallel “milestone” recognition method can be implemented using a field-programmable 
logic array (FPLA) since it has n inputs, n+1 outputs and w products. This is seen in 


the schematic below [Pe3]: 
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Figure 3. Serial-parallel PRBS to Natural Code Conversion (From: [Pe3]) 


The outputs B(n), B(n-1),...,B(1) yield the natural binary code for m*t which is 
the position of the recognized milestone Q(m). The output MS signals the control logic 
that a milestone as been detected. The control logic then stops the sequence of “back 
shifts” in the sequential code conversion method. The AGV stops traveling and 


computes the natural binary code of its position p by adding the binary representations of 


rand m*t using an n-bit adder [Pe2]. 
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For example, let n=5. Then the PBRS track is 2° -1=31 bits long. If the period is t =8, 
then the number of milestones is w= | (2° ia) 8 | =4. The Figure below illustrates this 


clearly [Pe2]. 


PSEUDORANDOM CODE TRACE 





Position 
O12 34 5 6 7 8 9 10 1112 13 14 1516 1718 1920 21 22 23 24 2526 27 28 29 30 Index, p 


Q0) am QQ) Q@) Milestones 


Figure 4. | PBRS Track with Milestones (From: [Pe2]) 


Each Q(i)(O <i <3) (boxed in with dashed lines) is a milestone and is associated with 


the following pseudorandom / natural code conversion scheme: 











| Milestone | Natural code position | Code decimal equivalent | 
00001 00000 0 
11101 01000 8 
O1111 10000 16 
11010 11000 24 

















Table 9. | Milestone-to-Natural Code Conversion (From [Pe2]) 


Assume the AGV’s initial position is 10011 (boxed in with solid lines). The 
AGV begins to travel one bit at a time towards position p = 0, comparing each n-tuple in 
parallel with all milestones. Once it reaches the milestone 01111, it obtains the position 
of this milestone to be m*t =16, where the natural binary code for this is 10000 (as seen 


in the Table above). Since the AGV needed to travel r =5 = 00101, bits to reach this 
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milestone, the absolute position of this n-tuple p =m*t+r=10000, + 00101, =10101, = 


21 bits which corresponds to the location of the n-tuple in the PBRS above. 


C. PERFORMANCE COSTS 


The performance cost of the previous example can be estimated as follows: 


Hardware cost: 4 words x 5 bits 


Time cost: 7 clock periods 


This cost is efficient when compared to a strictly parallel approach (31 words x 5 bits) or 
a strictly sequential approach (31 clock periods) [Pe2]. This comparison is demonstrated 
graphically below [Pe1]: 
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Figure 5. —_ Relative Time Performance Of Different Pseudorandom / Natural Code 
Conversion Methods (From [Pel1]) 


The following are equipment and temporal cost equations for the serial-parallel method 
employed by Petriu [Pel]: 
Equipment Cost = k,enumber of milestones + k, 


=k, (2"')/t | +k, a 
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Temporal Cost =k, +k,¢t (2.3) 
Total Cost = Equipment Cost + Temporal Cost = k, | (2"") /t] +k, +k,+k,et (2.4) 


where 


k, = equipment cost associated with each milestone 

k, = basal equipment cost for the serial back shift operations 

k, = basal temporal cost for a fully parallel solution 

k,= temporal cost associated with each back shift operation 

and k,, k,, k, and k, are constant for a given n and are functions of the technology 


being used [Pel]. The following graph illustrates this relationship: 
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Figure 6. — Serial-Parallel Code Conversion Costs as a Function of Distance (From [Pe1]) 


The total cost has a minimum at 7,,,and can be found by using the following formula: 


boyy = [« Toca (2.5) 


the optimal distance, t,,,, between milestones depending on the values of k, and k,. 


opt ? 
More parallelism is needed as the temporal cost begins to exceed the hardware costs 
k, >k,. As the measuring resolution increases, we will need more milestones to maintain 


the same code conversion speed [Pel]. 
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Iii. GREEDY DEBRUIJN SEQUENCE (GDBS) APPROACH AND 
RESULTS 


A. BACKGROUND 

In Petriu’s scheme, the AGV monitors its position on a binary track that is 
physically laid out on the floor. This track is modeled after a Pseudo-Random Binary 
Sequence (PRBS) that is generated using a linear shift register. Our scheme also uses a 
binary track laid out on the floor that is encoded using the greedy DeBruijn sequence. 
This sequence can be constructed from the concatenation of lexicographically ordered 
Lyndon words generated by the Necklace algorithm. The objective of this chapter is to 
demonstrate an alternate scheme of measuring the absolute position of any Automated 
Guided Vehicle (AGV) whose position is given by an n-tuple on a greedy DeBruijn 
sequence. 

There are two ways to measure the absolute position of an n-tuple depending on 
whether the position is above or below a threshold. The first method applies a direct 
computation for those n-tuples which are Lyndon words or which lie between two 
adjacent necklaces that are equal to or greater than the threshold. Since there are no 
missing even decimals above the threshold, we can directly calculate the position of any 
n-tuple above the threshold without any need of signposts or without running the 
Necklace Algorithm. The second method uses lexicographically ordered signposts that 
serve a similar function to Petriu’s milestones. These signposts are evenly distributed 
throughout all the necklaces below the threshold. The signpost data are 2n bits long 
where the first n bits correspond to the necklace designated as a signpost and the last n 
bits indicate the binary equivalent of the numerical value of the distance from a point of 
origin. We define the LSB of the greedy DeBruijn sequence whose position number is 
designated as “1”’as the LSB (or the rightmost bit) of the 000...0 n-tuple. 

B. ABSOLUTE POSITION MEASUREMENT ABOVE THRESHOLD 

In our approach, an AGV travels either to the left or right tracking at most n bits 
until the AGV’s position matches that of a necklace. We call a nearby necklace the 
position of reference (POR). While the AGV moves towards the POR it tallies the 


number of bits it has traveled. The subroutine below, called testforNecklace, is 
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embedded in another routine that searches the vicinity of the original n-tuple and tests 
each new n-tuple as the AGV moves to see if it is a necklace. The routine stops 
executing when testforNecklace returns a true Boolean value indicating a necklace was 


found. 


bool Necklaces::testforNecklace(long long bininput); 
{ 
bool foundnecklace = false; 
long long ndecimal = bininput; // binary decimal associated with n-tuple 
nstring = decToBin(ndecimal); // subroutine that converts decimal value into a binary 
// sequence 
shiftSequenceToNecklace(); // subroutine that shifts a binary sequence into its necklace 
// representative 
shiftdecimal = binToDec(); // subroutine that converts a binary string into its decimal 
// equivalent 
if(shiftdecimal = = ndecimal) // test to see if current n-tuple is a necklace 
{ 
foundnecklace = true; 
} 
else 
{ 
foundnecklace = false; 
} 
return foundnecklace; 


// end testforNecklace 


The key subroutine in testforNecklace is shiftSequenceToNecklace. Further details of the 


coding implementation can be found in Appendix B. 


The absolute distance of any n-tuple above the threshold is given by: 
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2" _n k-1 
p=2" -ne| 5 U] Sher pysb, (3.1) 
i=l 


where D is the distance of the most significant (leftmost) bit of the n-tuple from the 


point of origin (including the point of origin), n, is the decimal value of the POR, p,is 


the period of the i necklace generated by the Necklace Algorithm and b is the number of 
bits traveled to reach the necklace (where —b represents the AGV moving b bits to the 


left to reach a POR). For example, given the greedy DeBruijn sequence for n = 7 


14111111041111100%1111010411110001110110%111010041 1100101110000 1 101 100% 
1101010%1101000%1 10010041 10001041 1000004 10101004 1010000%1001000"1000000°0 


we determine the number of leading ones the threshold (highlighted above) to be 


p= a = H =3. Assume 1011101 (underlined above) is the AGV’s initial position 


(shown using an underline). According to the above subroutine, the decimal value of 
1011101 is 93 when the point of origin is the LSB position of the 0000000 n-tuple in the 
above sequence. Shifting the string to its necklace representative 

1011101 > 1110110 
and determining its decimal value we have 118 #93. So the AGV keeps moving to the 
left until the decimal values match. After moving five bits to the left it reaches the 
necklace 1110110, which is designated as the POR. Using the formula for D, we have 
the position of the 7-tuple 1011101 to be: 
ra, a8). Sa-n-s 





p=2"-14| 


D=128-—35+6+0+0+0+0-5=94 


which corresponds to the location of the sequence above. 


In the above example, we used n = 7, a prime number. When n is composite (e.g 


n= 6) we can see more clearly the effects of shortened periods. (See Table 10) 
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Necklace | Lyndon Word | Period | Decimal Equivalent 
111111 1 1 63 
111110 111110 6 62 
111100 111100 6 60 
111010 111010 6 58 
111000 111000 6 56 
110110 110 3 54 
110100 110100 6 a2 
110010 110010 6 50 
110000 110000 6 48 
101010 10 2 42 
101000 101000 6 40 
100100 100 3 36 
100000 100000 6 32 
000000 0 1 0 




















Table 10. | Necklaces and Lyndon Words for the Greedy DeBruijn Sequence 


We construct the greedy DeBruijn sequence given for n = 6: 


14111110%111100%111010%111000%110"1 101001 10010*110000°10*101000%100 
41000000 
6-1 


The threshold has p = al = a = | 2.5 | = 2 leading ones and is highlighted in our 





sequence. If the AGV’s current position is 100110 and we move three bits to the right to 
reach the POR 110010 we obtain 





6 1 
p=2-64|? 5] Denys 
i=l 


D=64—424+54343=33. 


This location is also confirmed by the above sequence. 
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It is not too difficult to determine the necklaces that have shortened periods for a 
given composite n before you run the Necklace Algorithm. First, determine the prime 
divisors of n. Partition the n-length sequence of all ones into subsequences whose 
lengths are the cofactors of these prime divisors (numbers resulting from the divisions by 
the prime divisors). We run the Necklace Algorithm on each of the partitions 
simultaneously until all of the necklaces have been generated for each of these identical 
partitions. Finally, we concatenate all identical partitions into a sequence of length n. 
This results in all the necklaces that have shortened periods. 

For example, let n =12. The prime divisors of 12 are 2 and 3. Dividing 12 by 2 
and 3 yields 6 and 4. Partitioning the 12-long sequence 111111111111 into the 
subsequences of lengths six and four yields 1111114111111 and 11114111141111. We 
then generate Table 11 for the cofactor d = 6 and Table 12 for the cofactor d = 4. 









































Necklace Period | Decimal 
1111114111111 > 111111111111 1 4095 
111110%111110 > 111110111110 6 4030 
111100%111100 > 111100111100 6 3900 
111010%111010 > 111010111010 6 3770 
111000%111000 > 111000111000 6 3640 
110110*110110 > 10110110110 3 1462 
110100110100 > 110100110100 6 3380 
110010*110010 > 110010110010 6 3250 
110000%110000 > 110000110000 6 3120 
101010101010 > 101010101010 2 2730 
101000101000 > 101000101000 6 2600 
100100100100 > 100100100100 3 2340 
100000* 100000 > 100000100000 6 2080 
000000*000000 > 000000000000 1 0 

















Table 11. Generating Shortened Necklaces for d = 6 
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Similarly, 























Necklace Period | Decimal 
11114111141111 > 111111111111 1 4095 
1110111041110 > 111011101110 4 3822 
1100%1100%1100 > 110011001100 4 3276 
1010*1010%1010 101010101010 2 2730 
1000*1000*1000 > 100010001000 4 2184 
000000000000 > 000000000000 1 0 

















Table 12. | Generating Shortened Necklaces for d = 4 


Note that since 6 and 4 have a greatest common divisor (g.c.d.) of 2, they share the same 
common factors of 1 and 2. Thus, they also share the necklaces whose periods are one 
and two, namely: 111111111111, 000000000000 and 101010101010. All the other 
necklaces for the two divisors are distinct. Since 3 is a divisor of 6 and not of 4, all 
necklaces with shortened periods of length 3 are contained in the list for d = 6 as well as 
d= 3. Since all necklaces are lexicographically ordered, we know which of the necklaces 
with shortened periods will precede the POR so we can make the necessary adjustments 
in calculating the AGV’s position. Of course, we can make things simpler by insisting 
that n be a prime number. 
C; SIGNPOST GENERATION BELOW THRESHOLD 

When measuring the absolute position of an n-tuple below the threshold we have 
amore difficult situation. Since there are “missing” even decimals to contend with as 
described in Chapter II we cannot apply a direct computation. We choose to select some 
of the necklaces below the threshold as “signposts” to contain embedded information 
about their distance from the origin. This procedure dominates absolute position 
measurement as n increases since the ratio of necklaces above the threshold to those 
below gets very small. This is confirmed through mathematical derivation and supported 
by experimental data illustrated in Table 16 and Figures 11 and 12. Although these 


signposts are distributed evenly within the DeBruijn sequence below the threshold, each 
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class will contain a different number of signposts since each class contains a different 


number of necklaces. We describe a placement of signposts for n = 7 in Table 13. 












































Necklace a Class | Period ea 
1111111 1 | 1 127 
1111110) 1111110] 6 7 126 
1111100)/1111100] 5 7 124 
1111010)/1111010} 4 7 122 
1111000;/1111000} 4 7 120 
1110110) 1110110] 3 7 118 
1110100) 1110100] 3 7 116 
1110010/1110010] 3 7 114 


















































1101100 1101100; 2 7 108 
1101010) 1101010] 2 7 106 
1101000 1101000; 2 7 104 
1100100)1100100| 2 7 100 
1100010,1100010) 2 7 98 
1100000); 1100000; 2 io 96 
1010100 1010100) 1 7. 84 
1010000) 1010000; 1 ei 80 
1001000 1001000) 1 vi 72 
1000000; 1000000] 1 7 64 
0000000 0 0 1 0 
Table 13. Necklaces and Signpost Insertion for n = 7 





Below the threshold (highlighted row), there are only 11 remaining necklaces within which 


we insert signposts. If we arbitrarily designate six necklaces as signposts (highlighted in 


31 


bold), we see class 2 has three signposts, class 1 has two signposts and class 0 has one 
signpost (the all Os n-tuple). Although we would not likely put the signposts in such close 
proximity when n is larger, this example serves to illustrate the point of equal placement. 
As n increases, the results are more dramatic and there are more necklaces between 
adjacent signposts. 

For any n, we generate the signposts as we run the Necklace Algorithm to 


generate the necklaces. To find the number of signposts, N,, we want to insert, we first 
need to calculate how many remaining necklaces, N,., there are below the threshold. We 


use the following formula: 








N,=@, [2 reo i} 
=Z 2" P41] (3.2) 
=Z,—N;, 


where p= Lal . Z, indicates the total number of necklaces for a given n and 


N, = 2" "+1 is the number of necklaces above and including the threshold. The 


N, 
N, -|*|, (3.3) 


where d < N,, d #0 and d represents the number of necklaces we want the necklace 


number of signposts, NV ,, is 





algorithm to generate between signposts on the list. We define the distance between 


signposts as N,. If nis prime, then NV, =d *n. Otherwise, since N, <d*n, we need to 


take into account necklaces that have periods less than n. In our example, 





p= 7 = = S ~3 and N, =Z,—2’*'-1=11 which agrees with the example above. 
2 


We choose d =1 yielding a total of six signposts and the distance between signposts is 
N, =1*7=7 bits. 
Recalling that we only concatenate the Lyndon words derived from the necklaces, 


we have the following greedy DeBruijn sequence for n=7: 
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1411111104111110041111010%1111000%1110110%1110100%1110010 “1110000 
411011001 101010*1101000%1 100100“1100010%1 1000001010100%1010000 
A1001000"1 00000040. 


The sequence is of length 2’ =128 with the signposts indicated in bold and the threshold 
highlighted. If we measure the point of origin from the zero bit at the LSB position of the 
DeBruijn sequence, then the number of bits needed to reach the MSB of our signpost is 
the absolute distance D. By beginning with a count of 2’ =128, we subtract the period 
of each necklace we generate to determine this distance, convert it into binary, and 


concatenate it to the necklace. This is illustrated in the Table below: 





Necklace | Distance (bits) Signpost 

1101100 | 71 = 1000111, ) 1101100%1000111 
1101000 | 57 =01110012 | 1101000°0111001 
1100010 | 43 =01010112 ) 11000100101011 
1010100 | 29 =00111012 | 1010100°0011101 
1001000 | 15 =0001111.2 | 1001000*0001111 
| 0000000 | 7 = 00001112 | 0000000°0000111 | 
































Table 14. Signposts for n = 7 


The signpost information has length 27, since it is a concatenation of the 
necklace designated as a signpost and the binary equivalent of its numerical distance 
from the point of origin. After we generate the signposts, the AGV can store them in its 
memory to use in determining its location. Like Petriu’s scheme, the AGV will find itself 
in an arbitrary position designated by a unique n-tuple. In Petriu’s implementation, the 
AGV must travel back towards the origin one bit at a time, using the window property of 
the PRBS to see if each n-tuple encountered matches a milestone. Since Petriu 
simultaneously uses a serial and a parallel operation to detect his milestones, this is an 


O(m) x O(/) operation, where m represents the distance between the n-tuple and the 
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milestone and / represents the number of milestones. Since the number of bits between 
milestones grows significantly as n increases, this can take a while and the work grows 
accordingly. 

Our signposts contain the same information as Petriu’s milestones. The main 
difference in the two schemes is that our signposts are lexicographically ordered, so their 
relative spatial relationship to each other is clear. The same cannot be said about Petriu’s 
milestones unless one performs the pseudorandom to natural code conversion described 
in Chapter III. 

Unlike in Petriu’s scheme, in our scheme each n-tuple examined is not compared 
with each signpost until a match is found. Instead, we search for a POR to compare 
against signposts so we can determine between which signposts the POR resides. Since 
the signposts are lexicographically ordered and organized by class number, we can 
narrow the search to a smaller subset of the list of signposts. As we compare the POR to 
our signposts, we keep track of the current and the last signpost we encounter on our list. 
When the first signpost is found whose decimal value is smaller than that of our POR, we 
stop the search. The AGV then runs the Necklace Algorithm from the last signpost 
whose decimal value is greater than the POR’s to the POR while simultaneously 
computing how many bits lie between them. This is accomplished by determining the 
period of each necklace generated and then summing the periods. This number is 
subtracted from the decimal equivalent of the distance information in the signpost. If the 
AGV needed to travel to the left to reach the POR, then we would subtract the number of 


bits traveled. Otherwise, we add them. We summarize this in the following formula: 
D=D,-), p,+b (3.4) 
i=j 
Dis the absolute distance of the n-tuple, Dis the distance of the i: signpost from the 


origin, p,is the period of a necklace between the j"" signpost and the m'" necklace that 


precedes the POR, and 5 is the number of bits traveled to reach the POR. 


For example, given the DeBruijn sequence for n=7, 
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1411111104111110041111010%1111000%1110110%1110100%111 001041110000 
A1101100%1 1010101101000%1 100100%1100010%1 100000*1010100%1010000*1001000 
4100000040 


If the AGV’s initial position is below the threshold, for example 0011001 (underlined 
above), we need to use the distance information contained in the signposts. If we choose 
to move to the left five bits, the AGV’s designated position will be the POR 1101000. 
Since this POR is already a signpost, we do need not run the Necklace Algorithm to find 
the separation distance between a signpost and the POR. Using the formula and the 


distance corresponding to 110100 from the Table above, 
D =57-—0-—5=52 bits from the point of origin 


There are no necklaces between the signpost and the POR (since they are the 
same) and the AGV needed to move five bits to the left. This matches the count we 
would obtain from the above sequence. If the POR is not a signpost, we need to search 
through all the signposts in the same class as our POR to find out between which 
signposts our POR resided. We then run the Necklace Algorithm to get the information 
needed to compute the absolute position. 

For example, assume we only choose to insert four signposts (indicated in bold) in 


our sequence below: 


1411111104111110041111010%1111000%1110110%1110100%111 00101110000 
4110110041 10101041 1010001100100" 1000101 100000“1010100%1010000 
41001000“10000000 


If we have the same POR as before, 1101000 is not a signpost. We compare this necklace 
to the list of signposts below 

1101100 

1100100 
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1010100 

1000000 
and we see that 1101100 < 1101000 < 1100100. Running the Necklace Algorithm from 
1101100 to 1101000, we find that they are 14 bits apart. Using formula 3.4 


m 


D=D,->.p,+b 


i=j 
D=71-(7+7)—-5=71-19=52 bits from the point of origin 
This corresponds with the previous result we had for this n-tuple. 

With more than one AGV on the DeBruijn track, we initially space them 
uniformly apart where their initial positions correspond to a signpost. Then, as each 
AGV begins to move we update their current locations by adding new signposts to the list 
or by changing the identities of the signposts. If the AGVs are not initially on signposts 
and not spaced evenly, we could move each of them to their respective PORs and easily 
find out how far apart they are from each other. Even without knowing the exact location 
of each POR, we can have a sense of where each AGV is on the DeBruijn track and 
where they are relative to each other since the PORs are lexicographically ordered. (A 
higher decimal value of the POR corresponds to a longer distance from the point of 
origin). This also leads to an optimized solution on the AGV placement since we can 
calculate (without moving them beyond the POR) which one is closest to a particular 
signpost or AGV. This is not possible with Petriu’s milestones, since in his formulation 
the milestones are not lexicographically ordered with respect to their distance from the 
point of origin. In addition, an AGV on Petriu’s track needs to move much further than n 
bits (when n is sufficiently large) to find its location to get a sense of how far it is from 
the point of origin. The problem becomes more expensive with more than one AGV on 
the track for a sufficiently large value of n. 

D. PERFORMANCE EVALUATION 

Petriu uses a pseudorandom-to-natural code conversion with his milestones. 

Depending on the value of n that is used to generate the PRBS of length 2” —1, the 


number of milestones can become large, or if we limit the number of milestones, the 
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distance between adjacent milestones can become large. Figures 7 and 8 indicate how 


the value of t., the distance between milestones and w, the number of milestones, 


opt ? 


grows with the value of n. 


Distance between milestones vs. n 
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Figure 7. Distance between Milestones vs. n 
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Number of milestones vs. n 
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Figure 8. | Number of Milestones vs. n 


Figures 9 and 10 compare the number of signposts versus n and the distance between 


signposts versus n. 
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Distance between singposts vs. n 
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Figure 9. _ Distance between Signposts vs. n 
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Figure 10. Number of Signposts vs. n 


Petriu’s milestones are fixed once they are chosen for a given n as they are 
designed into the hardware. Our signposts are stored in memory and, each time an AGV 
moves to a new location, we can determine a new POR and add a new signpost among 
existing signposts while the AGVs are in operation. We have the flexibility of adding or 
removing as many signposts as we want or changing which POR we choose to be a 
signpost. Table 15 compares the number of signposts and their distances between them 


with Petriu’s statistics. 
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Greedy DeBruijn Sequence Track 















































n 5 | 7 8 16 19 24 31 50 
2” | 32 | 128 | 256 | 65,536 | 524,288 | 16,777,216 | 2,147,483,648 | 1.12589990684262e+015 
Z, | 8 | 20 | 36 | 4116 | 27,596 699,252 69,273,668 22,517,998,808,028 
N, | 3 | 11 | 19 | 3859 | 27,083 695,155 69,240,899 22,517,965,253,595 
Ne 5 9 17 257 513 4097 32,769 33,554,433 

d 1 1 1 8 20 86 748 335,545 
D, 

N, 

Petriu’s Track 
Lot 
Ww 








Table 15. Statistics on GDBS vs. PRBS Performance 


D,= | (d —1)* n | refers to the distance between the LSB of one necklace to the MSB of 


t 
the adjacent necklace. The variable d = f= is chosen so that every d necklace could 
nN 


be selected as a signpost so that the distances between necklaces would match Petriu’s 


optimal distances. Both approaches (highlighted above) perform about the same. The 





: N.|. . 
number of signposts, J, -| 7 , 1s about the same as the number of milestones, w. 


We restrict our signpost selection to those necklaces below the threshold. However, the 


percentage of necklaces, %N,, above the threshold decreases rapidly as n becomes large. 


Given N, = Z,-2""'-1=Z —N,,,then 





gre 4+] 7 
Z Z Z 2" 
n 





4] 


n 


since p=|* |= ana ee Then as n> oo, 2" >>n aide > 1. 





n 


n 


Since VN, =Z,—N,, we have N, > 0as n—-> 00. This is also verified experimentally in 


Table 16 and in Figures 11 and 12. 





n 5 | 7 8 16 19 24 31 50 
YN, | 37.5 | 55 | 52.78 | 93.76 | 98.14 | 99.414 | 99.95 | 99.99985 


YN, | 62.5 | 45 | 47.22 | 6.24 | 1.86 | 0.586 | 0.05 | 0.00015 



































Table 16. %N,and %Nr vs. n 


Fraction of necklaces above threshold vs. n 


— %Nt 


— quadratic fit 
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Figure ll. %N,vs.n 
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Fraction of Necklaces below threshold vs. n 


— %Nr 


— quadratic fit 
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Figure 12. %N, vs.n 


This explains why the performance between Petiu’s approach and ours are 
similar. Our biggest savings occur during the AGV’s movement to its POR to test 
whether or not each n-tuple examined is a necklace. The number of comparisons is at 
most n, while in Petriu’s case, for n = 50, there is a potential of 67,108,665 comparisons. 
As mentioned previously, the AGV on Petriu’s track needs to travel at most 16,777,216 
bits in order to know its location. The AGV on our DeBruijn track would need to travel 
at most 50 bits. Of course, if the AGV does not find a signpost at the POR, then more 
computation is needed. If our POR lies between two signposts, for n = 50 the number of 
necklaces between adjacent signposts is d = 335,545 necklaces. Using the Necklace 
Algorithm to generate this number of necklaces for n = 50 is still less work than 
generating necklaces for n = 30 (since the number of total necklaces for n = 30 is 


364,724). Once we have our POR, we know in which class to begin our search. The 
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number of comparisons needed for our necklace will be less than 67,108,631 since all of 
50-1 

the signposts are spread out among p—1= 3 —1=23 classes. Assuming, as a 

baseline measure, that these signposts were evenly distributed among these classes we 


we) 
2 


would have =2,917,767 signposts over which to compare our necklace. In 


reality, some classes have more signposts than others. In Petriu’s scheme, comparisons 
against the milestones occur numerous times while in our scheme comparisons against 
our signposts happens only once. 
E. SHIFTING ONES SIGNPOSTS 

A way to measure the effectiveness of the signpost scheme is to measure the 
average distance between the signposts and where along the greedy DeBruijn track they 
are distributed. This gives an idea of how much computation the AGV will need to 
perform in order to locate its POR within a list of signposts. An alternate technique to 
generate signposts is termed the shifting ones signposts. The idea is to see if a different 
distribution of signposts throughout the DeBruijn track below the threshold could 
produce better results than evenly distributing them. We compare the average distance 
between the shifting ones signposts with the distance between evenly spaced signposts. 

First, we select the last necklace associated with a given class number. For 
example, with n = 5, 11000 is the last necklace associated with class 2. We then shift the 
“1” associated with the largest index one space to the right, each time checking that we 
still have a necklace. These necklaces will be associated with our signposts. For 


example: 


11000 
10100 
10010 Not a necklace 
10000 


For larger n, there are many necklaces that demonstrate this pattern but we only 


select those necklaces that are associated with a clump (a group of consecutive even 
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decimals that are “missing” between two adjacent necklaces). The reason for this is 
explained more clearly in Chapter IV. In the example above, if we select only those 
necklaces associated with “missing” decimal equivalents to be signposts, then the final 


list of signpost values 1s: 


11000 
10100 
10000. 


The value n = 5 is too small to make a big difference. However, for n = 15 when we run 
the Necklace Algorithm we obtain the following nine signposts out of a potential 97 


necklaces which demonstrate this shifting ones pattern: 


110000000000100: clump size = 1 even decimal “missing” 
101000000000000: clump size = 1 even decimal “missing” 
100100000000000: clump size = 3 even decimals “missing” 
100010000000000: clump size = 7 even decimals “missing” 
100001000000000: clump size = 7 even decimals “missing” 
100000100000000: clump size = 127 even decimals “missing” 
100000010000000: clump size = 63 even decimals “missing” 
100000000000000: clump size = 63 even decimals “missing” 


The clump size refers to the number of “missing” decimal equivalents. These signposts 


are not evenly spaced apart as can be seen from the Table 17: 

















Necklaces Distances Between Necklaces 
31359 
110000000000100 
1085 
101000000000000 
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Necklaces Distances Between Necklaces 

183 
100100000000000 

60 
100010000000000 

20 
100001000000000 

15 
100000100000000 

15 
100000010000000 

15 
100000000000000 

15 
000000000000000 














Table 17. Distances between Shifting Ones Signposts (weight = 1) for n = 15 


The average distance between signposts is 3640 bits and the first distance corresponds to 
the distance from the MSB of the DeBruijn sequence to the MSB of the first signpost. 
The distance between our evenly spaced signposts is 90 bits and we need 295 of them. 
The DeBruijn sequence itself is of length 2'° = 32,768 bits. The beginning gap is 31,359 
bits long so the percentage of the DeBruijn sequence covered by signposts is 


ELS — 31,359 


*100 =4.3% 
32,768 


which is very inefficient since only a small percentage of the DeBruijn track is covered. 
If we decrease our searching range by increasing the weight to 2, we obtain 39 


(containing a weight of one or two) signpost ts out of a potential 361 necklaces: 
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Signposts 


Distances between Signposts 



















































































23,336 
111000000001100 
4778 
110100000000100 
1620 
110010000000100 
725 
110001000000100 
375 
110000100000100 
45 
110000011000000 
165 
110000010000100 
120 
110000001000100 
i 
110000000100100 
45 
110000000010100 
45 
110000000001010 
30 
110000000000100 
325 
101010000000000 
230 








101001000000000 
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Signposts 


Distances between Signposts 



















































































120 
101000100000000 
75 
101000010000000 
45 
101000001000000 
30 
101000000100000 
15 
101000000010000 
15 
101000000001000 
15 
101000000000100 
15 
101000000000000 
63 
100100100000000 
45 
100100010000000 
15 
100100001000000 
15 
100100000100000 
15 
100100000010000 
15 








100100000001000 
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Signposts Distances between Signposts 

15 
100100000000000 

15 
100010001000000 

15 
100010000100000 

15 
100010000010000 

15 
100010000000000 

15 
100001000010000 

5 
100001000000000 

15 
100000100000000 

15 
100000010000000 

15 
100000000000000 

15 
000000000000000 














Table 18. Distances between Shifting Ones Signposts (weight = 2) for n = 15 


The average distance is 840 bits. The DeBruijn sequence covered by signposts is 


— — 23,336 


*100 = 28.78% , 
32,768 


which is a significant improvement. However, the percentage of the 
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DeBruijn sequence covered by equally spaced signposts is 
ny 
[Poe sro. (3.5) 


For n = 15, we have 


2'° -15*129 
— 


; }100=94.10%, 


which is much more efficient. The following Table summarizes the performance of the 


shifting ones signposts as the weight is increased: 


























| Weight | Distance to first signpost (bits) | Number of Signposts | % of DB sequence | 

1 31,359 9 4.3 

2 23,336 39 28.78 
3 13,701 116 58.19 
nl 6976 244 78.71 
5 2866 370 91.25 
6 1921 432 94.14 
i 1921 447 94.14 
8 1921 450 94.14 














Table 19. | Number of Shifting Ones Signposts vs. % of DeBruijn Sequence for Various 
Weights 


The gap between the first signpost and the beginning of the DeBruijn sequence 
levels off to 1921 bits yieldig the highest efficiency of 94% with at least 432 shifting ones 
signposts. Evenly spacing the signposts throughout the DeBruijn sequence achieves the 
same efficiency but with 295 signposts. Clearly, the evenly spaced signpost scheme has a 


better performance. 
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IV. CONSIDERATIONS AND FUTURE WORK 


A. BACKGROUND 
In Petrui’s scheme, the AGV needs to examine every n-tuple sequentially until it 
reaches a milestone. There is no pattern in the sequence it traverses. In our scheme, the 
AGV needs to initially examine at most n bits sequentially, and then it can examine n bits 
simultaneously as it compares the POR to the signposts because of the pattern inherent in 
the list of necklaces. Although our scheme is efficient in finding the absolute position of 
an n-tuple, we have attempted to improve on its efficiency by a couple of methods. The 
first method involves trying to find a pattern in the number of missing necklaces that 
allows us to calculate the number of “missings” in a given class for any n. The second 
method involves trying to map the necklaces of length n to a class of a larger group of 
necklaces of length n+k, where k is some class number for this larger group of necklaces. 
Rather than scan through the list of signposts one at a time to find the location of a POR, 
either method would allow us to make calculable leaps over a group of signposts and 
speed up the search problem. Determining a method of computing the position of an n- 
tuple in a greedy DeBruijn sequence would eliminate the need to use the Necklace 
Algorithm to calculate the position. It would also eliminate the need for having a 
physically laid out binary track for the AGV to follow. The AGV would know the 
location of the n-tuple that designates its position without having to travel to its POR on 
the binary track. Although there are readily apparent patterns, finding a mathematical 
relationship to describe these patterns has been challenging and would be an appropriate 
topic for future work. 
B. RECURRENCE RELATION FOR MISSINGS 
h 


A linear recurrence relation is of the form h,,, =c,_,h,,,_,+...+¢,h,,, where the 


n-1 i+] 

c,'s (l<i<n-—l) are integers and the h,'s (1<i <n) (in our presentation) are positive 
integers representing the number of “missings” for a given group of necklaces of length 
i. We attempt to predict the number of “missings” for a given value of n knowing the 


number of “missings” for smaller values of n. 
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We ran the Necklace Algorithm and collected and organized the following 
information on the number of “missings” per class for a given n. (The class number is 


determined by k = n-q). 








































































































q_ |\n=2\n=4|\n=6|n=8 \n=10 n=12 |n=14| n=16 |n=18 | n=20 | n=22 
0;0};0;/0;0/ 0); 0 | 0 0 0 0 0 
1;0;0;0;/0;] 0/0 | 0 0 0 0 0 
2/0/0/0]0}] 0 | 0] 0 0 0 0 0 

3 0;O0;0/;/ 0; 0 | 0 0 0 0 0 

4 3|;0;0;] 0/0 | 0 0 0 0 0 

5 3|0]0/]0 /] 0 0 0 0 0 

6 15|5 | 0 | 0 | O 0 0 0 0 

7 25| 5 | 0 | 0 0 0 0 0 

8 63| 31 | 5 | 0 0 0 0 0 

9 114} 32 | 5 0 0 0 0 
10 255 | 155 | 33_| 5 0 0 0 
11 482 | 174 | 33 5 0 0 
12 1023] 719 | 180 | 33 5 0 
13 1985| 846 | 181 | 33 5 
14 4095| 3156 | 897 | 182 33 
15 8050 | 3911 | 916 182 
16 16383)13469) 4256 | 922 
17 32440} 17501 | 4394 
18 65535) 56458 | 19531 
19 130307} 76561 
20 262143) 233726 
21 522478 
22 1048580 





Table 20. | Number of “Missings” for Even n 


Notice that as n increases the number of “missings” per class approaches a steady state 
value. From the Table above, our steady state values (highlighted in bold) are 5, 33 and 
182. 
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q |n=1n=3 |n=5\n=7 n=9 n= N=13 n=15| n=17 | n=19 | n=21 
0;o0;0;0;0;0/; 0 |] 0 | 0 0 0 0 
1;0/;0;0;0]0] 0/0 /| 0 0 0 0 

2 0O;O;/0/0}]0|0 0 0 0 0 

3 1/0/;0/0}] 0 |0 0 0 0 0 

4 2,0;0/ 0 | 0 | 0 0 0 0 

5 d [2 | 0}; 0} 0) 0 0 0 0 

6 12}2]{0 /|0 | 0 0 0 0 

7 31/13] 2 | 0 | 0 0 0 0 

8 55| 14 | 2 | 0 0 0 0 

9 127| 72 | 14 | 2 0 0 0 
10 238 | 78 | 14 | 2 0 0 
11 511 | 340 | 79 | 14 2 0 
12 984 | 389 | 80 14 2 
13 2047/1515] 408 | 80 14 
14 4003| 1834| 414 80 
15 8191} 6546 | 1970 | 415 
16 || 16174} 8312 | 2021 
17 32767| 27642 | 9158 
18 65044 | 36708 
19 131071115037 
20 260975 
21 524287 


Table 21. | Number of “Missings” for Odd n 


The steady state values in the Table for n odd are 2, 14 and 80. There is an additional 
pattern that is shared by both Tables. These values are verified by listing the necklaces 
for a given n and a given class k = n-q. For example, if n = 7 we have the following 


information in Table 22: 


a 





















































Necklace Class Decimal Equivalent 
1111111 7 127 
1111110 6 126 
1111100 =) 124 
1111010 4 122 
1111000 4 120 
1110110 3 118 
1110100 3 116 
1110010 3 114 
1110000 3 112 
1101100 

1101010 

1101000 

1100100 

1100010 

1100000 

1010100 1 84 
1010000 1 80 
1001000 1 72 
1000000 1 64 
0000000 0 0 

















Table 22. Necklaces and Their Classes for n = 7 Used in Analyzing ““Missings” 


Note for class 7-5 = 2 there are two missing decimal equivalents: 110 and 102, 
corresponding to the strings 1101110 and 1100110, thus the entry (highlighted) in the 
Table above of 2 “missing”. When considering class 9-6 = 3 (for n = 9) and class 11-7 = 


4 (for n =11) there are also 2 “missings”: 111011110 and 111001110 (forn =9), 
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11110111110 and 11110011110 (for n = 11). With knowledge of how many “missings” 
there are between necklaces, we can use the following formula in determining how many 


necklaces, NV, there are between any two necklaces 
n, -n 
n,-{ - )-yg-1 (4.1) 


where n, ,n,, are the decimal value of the necklaces and n, 2n, , and N,, represents the 
number of “missings” between these two necklaces (N,, =0 if n, =n,,). For example, 


given the necklaces 1110000 and 1100000, we have 


N= (ee) —2-1=5 necklaces , 


n 


which corresponds with the Table above. Without knowledge of the number of 
“missings”, we would have needed to run the Necklace Algorithm in order to determine 
this information. 

Although the steady state values differ whether n is even or odd, they do share an 
additional common pattern. Comparing the first few diagonal elements as we approach 


the steady state values we have: 


For n even: 

3, 3,5, 5,5, ... steady state value = 5 

15, 25, 31, 32, 33, 33, 33, ... steady state value = 33 

63, 114, 155, 174, 180, 181, 182, 182, 182, ... steady state value = 182 


For n odd: 
1, 2, 2, 2, ... steady state value = 2 
7, 12, 13, 14, 14, 14, ... steady state value = 14 


31, 55, 72, 78, 79, 80, 80, 80, ... steady state value = 80 


If we form the difference sequence between adjacent values we have: 


aD 


For n even: 

2, 0, 0, 0, ... 

10, 6, 1, 1, 0, 0, ... 

51, 41, 19, 6, 1, 1, 0, 0, ... 


For n oda: 

1, 0, 0, 0, ... 

5; 1, 1, 0,0; Opa 

24, 17, 6, 1, 1, 0, 0, 0, ... 


Based on the difference sequences, we can conjecture what the steady state value should 
be for a given n, class number and the last few values approaching the steady state. In the 
case of n = 24 and class number = 24 —17 = 7, the steady state value should be 922+2 = 
924. For n = 23 and class number = 23 — 16 = 7, the steady state value should be 415 + 1 
= 416. 

Running the data for higher values of n, we notice the difference sequences for n 
odd or even are the same, namely: 1, 1, 6, 19, 51, 138, ... Comparing these values to an 
online integer sequence database [Sl] to see if this pattern matched any other 
mathematical structure, we found there was no match. Evidently, this problem has not 
been studied previously. More work needs to be done to uncover the meaning of this 
pattern since it unifies the data for both n odd and even. 

C, MAPPING OF SUBSEQUENCES OF NECKLACES 

A second approach to understand the missing n-tuples involves trying to 
enumerate the number of necklaces in a given class k = n-q, by creating a mapping of 
these necklaces to those of length n-k, where the total number of necklaces was known. 
If we know a lot of information (such as number and location of “missings’’) for 
necklaces of length n-k, we want to know how much information we can predict about 
necklaces of length n, class k. For group of missings associated with steady state, there 


seems to be a lot of predictability. 
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Define an “m-missing” to be a binary string in the group of “missings “associated 
with steady state value m for a particular n. For example, the “2-missings” for n = 7 andn 


= 9 is shown in Table 23. 





n =7 “2-missing” (Class 2) | n = 9 “‘2-missing” (Class 3) 
11011102 = 110j0 1110111102 = 47810 
11001102= 10210 1110011102 = 46240 

















Table 23. Comparison of n = 7 and n = 9 “2-missings” 


The binary and deciemal values are shown. Note that the “2-missings” for n = 9 is 
obtained by adding a | to the two longest runs of 1s for n = 7 (in this case there are only 


two runs in which to add a 1). 




















n =8 “5-missing” (Class 2) | n = 10 “‘5-missing”’ (Class 3) 
11011110) = 22246 11101111102 = 95819 
110111002 = 220j0 11101111002 = 95640 
110101105 = 21446 11101011102 = 942; 
110011102 = 20610 11100111102 = 92619 
110001105 = 198j6 11100011102 = 91040 














Table 24. | Comparison of n = 8 and n = 10 “5-missings” 


In both Tables the length n increased by two and | bit was added to each of the longest 
run of ones. Since we have a one-to-one mapping between both sets of “missings”’, the 
number of “missings” does not change and that is why we have steady state values for the 


“missings”. 


There is another means of evaluating a mapping from a set of missings for length n-k to 
another set of missings of length n. Consider the clump that is generated as the Necklace 


Algorithm transitions from one class to another while generating necklaces. The last 


my 


element of a class k will have k ones followed by all zeros: 111...10....0. The next 


kones 


necklace produced will be of the form 111...10111...10...0111...10 11..0 where ris the 


k-lones k-lones k-lones n—r(2k-1) 





number of times 111...10 occurs in the sequence. Using the following formula, we can 


k-lones 
find the number of missings between two consecutive necklaces generated by the 


Necklace Algorithm: 





N {ls —Na, 


=e 4.2 
5 (4.2) 


N,, is the number of missings and NV, , N, are the decimal values of the necklaces. We 


can then determine the clump size, or number of necklaces, between the smallest 
necklace of class k and the laragest necklace of class k-1. If we remove the first k bits 


from the left of the necklace111...10111...10...0111...10 11..0 , then we obtain a smaller 


k—lones k—lones k—lones n-r(2k-1) 





version of this necklace that is derived from the necklace 111...10....0 of length n-k. 


kones 
Knowing the size of the first clump for the n-k case should give us insight into knowing 


the size of the first clump for the n case. For example, let n = 32. Then for k = 6 we have 
11111100000000000000000000000000. 

Applying the @-Algorithm, we obtain 
11111011111011111011111011111010. 

Removing the first 6 bits from 11111011111011111011111011111010 we obtain 


11111011111011111011111010 (m7 =26). 
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This is derived from (using the Necklace Algorithm) 
11111100000000000000000000 
The following Table summarizes the number of missings between these consecutive 


necklaces for n = 32, 26, 20, 14 and 8. 














n= 32 n= 26 n=20 n=14 n=8 
61’s 
followed by 4227858432 66060288 1032192 16128 252 
all Os 
g! 4226793210 66043642 1031930 16122 250 
No. 532610 8322 130 2 0 
Missings 




















Table 25. | Comparison of number of “missings” for n = 32, 26, 20, 14 and 8. 


If we take the following ratios 





oe Ey oe ol eee ee oe ea eo 
8322 130 ze 


it seems as n increases by 6, the number of missings for the first clump increases by a 
factor of 2°. This type of information can be useful in reducing the number of 
computations needed when generating necklaces. More work needs to be done to 
understand the nature of this mapping between various classes for necklaces of a given 
length n and other families of necklaces of length <n. 
D. CONCLUSION 

Improvement in Petriu’s scheme involved improving the hardware 
synchronization design. Improving our scheme is of a mathematical nature. Gaining 
more theoretical understanding of these issues will allow us to make significant 


computational gains. This potential benefit is not possible through Petriu’s scheme. 
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APPENDIX A: PSEUDO-RANDOM BINARY SEQUENCES 


A. SHIFT REGISTER GENERATION OF A PSEUDO RANDOM BINARY 
SEQUENCE (PBRS) 


A Pseudo-Random Binary Sequence (PBRS) is a shortened De Bruijn sequence 
since it contains all possible 2" n-tuples except the all zeros sequence [Mi]. So the length 
of the PBRS is 2"—1. A PBRS contains 2”" ones and 2"'—1 zeros. There are 2”' odd 
numbers (binary numbers ending in 1) and 2” ' —1even numbers (binary numbers ending 


in 0) between 1 and 2”". So the total is 
271 42"'-1= 52" +2”")—-1=2"~—1 total n-tuples 


The number of runs (consecutive sequence of identical integers) of ones and zeros are 
approximately the same. For any given n, we have % of the runs have length 1, % of the 
runs have length 2, 1/8 have length 3 and 1/16 have length 4, etc. as long as the fraction 


makes sense [Go]. 


One can generate a PBRS using a linear shift register that is modeled using a primitive 
polynomial h(x) of degree n [Ma]. A primitive polynomial is one that is irreducible and 
has maximum period. The coefficients of the primitive polynomial come from the field 


of Z, = {0,1}. An example is the following: 


h(x) =1x* +0x° +0x? +1x' 4+1x° = xt +x 41 (A.1) 


The diagram below gives an illustration of a linear shift register that generates a maximal 
length sequence of length 2*—1. Each stage in the register can contain a value of “O” or 
“1”. The only non-zero coefficients of h(x) are those corresponding to x*, x' and x°. The 
values in the stages corresponding to x' and x° are the only ones that are “tapped” to be 
added. Since x* = x' + x°, the value associated with x*is inserted into the stage associated 
with x°at the next time unit when all of the constants shift to the right. The all zeros 


input is excluded so that we do not produce a continuous sequence of zeros. Starting 


with an initial state of 1 0 0 0 loaded into the shift register in Figure 13, we generate the 
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various state values as a result of adding the stages corresponding to x' and x”, and then 


shifting the contents of the stages one step to the right [Ma]. (See Table 25). 


OUTPUT 


0*x? 0*X? 1*x! 1*x? 





Figure 13. Feedback Shift Register Corresponding to x +x +1 (From [Ma]) 










































































State 
State Number | MSB LSB Output 
x x e x” 

0 1 0 0 0 0 
1 0 1 0 0 0 
Z 0 0 1 0 0 
3 1 0 0 1 1 
4 1 1 0 0 0 
5 0 1 1 0 0 
6 1 0 1 1 1 
Z 0 1 0 1 1 
8 1 0 1 0 0 
9 1 1 0 1 1 
10 1 1 1 0 0 
11 1 1 1 1 1 
12 0 1 1 1 1 
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State 
































State Number | MSB LSB Output 
a a a: 5 ae 
13 0 0 1 1 1 
14 0 0 0 1 1 
15 1 0 0 0 0 
16 0 1 0 0 0 


























Table 26. 16 Feedback Shift Register States Corresponding to ae al (From [Ma]) 


Note that there are 15 states corresponding to all the possible 2*—1 different non-zero 
binary 4-tuples. Since states 15 and 16 are just repeats of states 0 and 1, we have 
generated a cycle of length 2*—1. The value for Most Significant Bit (MSB) column in a 
row is just the sum of the LSB value (the value for x°) and the x' value from the 
previous row. The output corresponds to the last column of states (the Least Significant 
Bit (LSB) position). Since the shift register essentially generates a cycle of length 2” —1, 
any one of the states we load into the shift register will produce a cyclically shifted 


version of the output of various initial states as can be seen in Table 26 [Ma]. 






































STATE 

S112 OUTPUT 

1 0 0 90 O000100110101111 
0 1 0 0 001001101011110 
0 O 1 0 010011010111100 
1 0 O 1 100110101111000 
1 1 0 0 001101011110001 
0 1 1 0O 011010111100010 
1 oO 1 1 110101111000100 
0 1 0 1 101011110001001 
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STATE 
OUTPUT 

3 2 1 0 

KN AE ee 

1 O 1 O 010111100010011 
1 1 0O 1 101111000100110 
1 1 1 O 011110001001101 
1 1 1 1 111100010011010 
0 1 1 1 111000100110101 
0 1 1 110001001101011 
0 0O O 1 100010011010111 














Table 27. 16 Output Sequences Corresponding to x* + x + 1 (From []) 


Thus there are 2*—1 different pseudo-random binary sequences, each one corresponding 
to a state of the linear shift register. (We can generate a sequence of length 2”, called a 
De Bruijn sequence only if we use a non-linear shift register [Go]). The important point 


to make is that the sequence is a result of how we define h(x). If we change h(x) to be 


another irreducible polynomial x* + x° +1, then we reverse the sequences in the previous 
Table. Appropriate polynomials exist for every value of n [Ma]. 
B. FINITE FIELDS AND SHIFT REGISTERS 
Let a, b and c be elements of a set F with two binary operations: addition and 
multiplication. Then F is a field if the following properties hold [Ma]: 
(i) a+b=b+a 
(ii) ab = ba 
(iii) at+(b+c)=(at+b)+c 
(iv) a(bc) = (ab)c 
(v) a(b+c)=ab+ac 
Properties (1) and (ii) are the commutative properties, properties (111) and (iv) are the 
associative properties and property (v) 1s the distributive property. Further, there must 


exist elements 0, 1 in the set and for every non-zero element (a # 0), the additive and 
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multiplicative inverses,—a anda'' respectively, must also exist. Then the following hold 
[Ma]: 

(vi) O+a=a 

(vii) (-a)+a=0 

(vill) Oa=0 

(ix) la=a 

(x) (a')a=1 

A finite field is one that contains a finite number of elements. A field of g elements is 
known as a Galois field and is denoted by GF (q) where gq is a power of a prime number. 
An example of a simple field is the integers modulo p , GF (p), where p is a prime 
number. The elements of this field are {0, 1,.... D—- 1\ and addition, subtraction, 
multiplication and division (by non-zero elements) are carried out modulo p . So if 


p=2, thenGF(2) = {0,1} =Z,. We can construct a field with p” elements, where n is 


any integer and p is a prime number, by selecting an irreducible polynomial h(x) of 
degree n. The elements of this field will all be polynomials in x of degree <n with 
coefficients from GF (p). We will illustrate with p=2 and n=4. We construct a field 
of 2* =16 elements. In this case, we let the irreducible polynomial be h(x) =x*+x +1. 


This leads to the results in Table 27 [Ma]. 




















Polynomial | Polynomial Coefficients | Power of ¢ 
1 0001 €o =] 
x 0010 é! 
= 0100 & 
x 1000 ee 
x+1 0011 é 
Xv +X 0110 e 
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Polynomial | Polynomial Coefficients | Power of ¢ 
+x? 1100 il 
x +x41 1011 e! 
x +1 0101 e 
x+x 1010 a 
x +x+1 0111 a 
4x +X 1110 rage 
4x7 txt] 1111 ee 
+x +1 1101 a 
x +1 1001 oe 

1 0001 &S =1=£° 

















Table 28. Polynomials Associated with 16 States of the Shift Register (From [Ma]) 


Note that the coefficients of the polynomials correspond to the different possible states of 
the shift register. (The 16" state, not appearing in the sequence, is the all zeros state). 
Note that the third column of the Table has corresponding powers of €. We denote é as 
the primitive element that is a root of the primitive polynomial h(x) = x*+x +1. If we 
have h(x) =x*+x +1=Oandé is aroot, then €*+&+1=0 or €* =&+1. From this 
relation we can derive all other entries in the Table. We can then find the powers of any 
polynomial or the products of any of the polynomials above simply by performing 


exponential addition on the powers of ¢€. 
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APPENDIX B: CODE 


A. BACKGROUND 

The most important subroutine is the runNecklaceAlgorithm code. The algorithm 
of runNecklace Algorithm had to be modified from Matty’s version [Mat]. Since the 
number of necklaces grows prohibitively large, as n gets large, it was necessary to gather 
statistical information on necklaces (such as number of “missings’’) and insert signposts 
while executing the Necklace Algorithm to generate necklaces. The level of detail 
analysis is also limited by the amount of memory needed to store information on 
“missings”, clump sizes and other pertinent information. The runNecklaceAlgorithm is 
also limited to generate necklaces of a given binary string length n at a particular time. 
So comparison of data for different binary string lengths can only occur by storing the 
data separately. Fortunately, it is not necessary to generate all the necklaces using the 
necklace algorithm so detailed analysis can be performed on a small range of necklaces. 
The Automated Guided Vehicle (AGV) determines its position in a greedy DeBruijn 


track by using this feature. 


B. HEADER FILES 
1. Necklace Header File 


J [PS 7 2 3 2 2g ae 2g ae ae ae ae ae he he fe he fe fe he 6 6 6 26 28 28 28 28 28 28 28 28 28 28 28 28 2 2 2 2 2 2 2 2 2 2 2 2 2 2 28 2 2 2 ee he he fe he fe fe fe fe he ie ie 2c 2s is 


// Necklaces.h 

// 

// LT John Ortiz 

// 

// Project: Thesis 

// Operating Environment: Windows XP Home 
// Compiler: Visual Studio .NET 

// Date: 

// Description: 

// 

PS 7 2 2 2 ag ae ae ae 2g ae ae he he he he he fe fe he he 26 26 216 28 28 28 2S 28 28 28 28 28 28 28 2 2 2 2 2 2 2 2A 2 2h 2 2 2 2 28 28 2 3 2 ee he fe fe he he fe fe he 2c 2k =e ie 2s is 
#ifndef NECKLACES_H 

#define NECKLACES_H 


67 


J [7B FR AH 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 he fe fe fe he he he Ie ie ie Fie ie Ie it 


// Class: Necklaces 
// Purpose: 


J [78 FR AR AR A AB A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 2 ee he he fe he he he he he Ie ie ie ie Ie is 


class Necklaces { 


public: 
const static int MAX COLUMNS = 400; 


// Keep in mind the limits of memory when displaying missings. 
// You might not display them all if there are too many. 
// Adjusting the threshold and percentage may help you see most of them. 


const static int DEFAULT_THRESHOLD = 0; 


const static int DEFAULT_PERCENT = 0; 
const static int MAX LENGTH = 40; 
const static int OUTPUT_CELL_SIZE = 4; //output size 


//**m = 14 is the largest dimension I can use here without running out of memory 


J [78 FR FR A 2 2s 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 he fe he fe he he he ie ie ie ie ie ie 


// Method Necklaces--Default Constructor 
// Return value none 

//Parameters none 

// Purpose 

J [7B FR AR A A A 2 2s 2 2 2 2 2 2 2 2g 2 2 2 3 2 2 ae ee he he he fe he fe ie ie ie ie ie 2k 


Necklaces(); 


J [7B FR A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee he fe he fe fe fe Ie ie ie ie 2c ie 2s 


// Method Necklaces-- Constructor 

// Return value none 

//Parameters int size, int thresh, float per 

// Purpose 

J [78 7 FR A A A 2 2 2 2g 2 2 2 2 2 2 2 2 2 2 2 2 3 2 ae ae he he fe he fe he fe ie ie 2k ie ie Ie is 


Necklaces(int size, int thresh, int per); 


J [FR FA i 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2h 2 2 2 ee he he fe fe fe he he ie ie Fie ie Ie it 


// Method initNecklace 

// Return value decimal equivalent 

//Parameters int ones 

// Purpose to create a sequence consisting of a certain number of leading ones 
W with the remaining sequence consisting of zeros 

J [78 FR A AR A A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 2 ee he he he he fe he fe Ie 2k iE 2k 28 Ie ik 


long long initNecklace(int ones); 
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J [7 FR A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2h 2 2 2 6 he fe he fe fe he he he ie ie ie ie 2Ie it 


// Method printSequence 

// Return value none 

//Parameters none 

// Purpose displays necklaces, theta generated and unobserved binary sequences 
J [7 FR A AR A A 2 2 2 2 2 2 He 2 2 2 2 2 2 2 2 2 3 2 ee he he fe he fe he fe ie Ie ie ie is Ie is 


void printSequence(); 


J [RFR A A 2 2s 2 2 2 2 2 2 He 2 2 2 2 2 2h 2 2h 2 2 2 2 2 fe fe fe he fe he he ie 2g ie ie Ie is 


// Method binToDec 

// Return value none 

//Parameters none 

// Purpose converts binary number of a sequence into its decimal equivalent 
J [78 FR A AR A A 2 2 2 2 2 2 He 2 2 2 2 2 2 2 2 2 3 2 ee he he he he fe he he ie 2c ie ie ie Ie is 


long long binToDec(); 


J [78 FR FR A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee fe he fe he fe ie ie ie ie ie ie 2s 


// Method decToBin 

// Return value none 

//Parameters int input 

// Purpose converts binary number of a sequence into its 

// decimal equivalent and returns a pointer to binary sequence 
J [78 FR AR A A A 2 2s 2 2 2 2 2 He 2 2 2 2 2 2 2 2 ee he he he fe fe fe he ie 2c ie 2k ie 2k 


void decToBin(long long input); 


J [78 FR AR A A 2 2 2s 2 2 2 2 He 2 2 2 2 2 2 3 2 3 ee he he he fe fe he fe ie ie ie ie ie ik 
// Method createUnobservedSequences 

// Return value none 

//Parameters int decl, int dec2 


// Purpose generates all even binary sequences between any two even sequences 
J [7 FR A AR A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 2 ee he fe fe fe fe he fe ie 2k ie ie 2s Ie ik 


long createUnobservedSequences(int a, int b); 


J [78 FR AR A A A 2 28 2 2 2 2 2 2 2 2 2 2g 2 2 2 2 2 ee he he fe fe he he ie 2k ie 2k 2s Ik 


// Method findDivisors 
// Return value none 
//Parameters int num 


// Purpose find all divisors for an integer n 
J [78 7 AR A A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 ae 2 he he he fe fe fe fe ie ie ik ie ie 2k 


void findDivisors(int num); 


J [78 A AB A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee he fe fe he fe he fe he Ie ie ie it Ie is 


// Method printDivisors 
// Return value none 
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//Parameters int num 
// Purpose prints all divisors for an integer n 
J [78 7 AR AB A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 2 ee he he fe he fe fe fe 2k 2k ie 2k ie Ie is 


void printDivisors(int num); 


J [RFR A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee he he fe fe fe he he ie ie ie ie 2S it 


// Method findGCD 

// Return value int gcd 

//Parameters int numl, num2 

// Purpose find greatest common divisor (gcd) for two integers 
J [7 FR A AR A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 3 2 fe fe fe fe he he fe he ie ie ie 2k Ie is 


int findGCD(int num1, int num2); 


J [78 FR A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 he fe he fe he fe he ie ie ie ie ie 2S 


// Method eulerTotient 

// Return value long long Zn 

//Parameters int num 

// Purpose Calculate number of necklaces for a given n 
J [78 7 AR A A A 28 2s 2 2 2 2 2 2 2 2 2 2 2 2 2 ae 3 ee he fe fe fe he fe ie ie ie 2k ie 2k 


long long eulerTotient(int num); 


J [FRR FF A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee he he he fe fe he he he ie ie ie ie 28 is 


// Method generatePower2 

// Return value long long power2val 
//Parameters long long exponent 

// Purpose generates 2 raised to any integer 

J [78 FR A AR A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ae ee he fe he fe he fe ik 2k ie ie ie Ie is 


long long generatePower2(long long exponent); 


J [7B FR A 2 2 2 2 2 2 2 2 2 2 2 2 2 2h 2 2 2 2 2 2 he fe he fe he fe Ie ie ie ie ie ie 2s 


// Method runNecklaceAlgorithm 

// Return value none 

//Parameters long long initial, long long final, long long input, bool numflag 

// Purpose executes necklace algorithm 

J [78 FR AR A A A 28 2 2 2 2 2 2 2 2 2g 2 2 2 3 2 3 ee fe fe fe fe fe he fe ie ie ie ie ie Ik 

void runNecklaceAlgorithm(long long initial, long long final, long long input, bool 
numflag); 


J [7B FR AR AB A A 2 2 2 2 2 2 2 He 2 2 2 2 2 2 2 2 2 2 2 2 fe fe he fe fe he fe ie Ie ie ie ie 2S is 
// Method inputBinaryString 

// Return value none 

//Parameters none 


// Purpose checks whether user input binary string is valid 
J [78 FR A RA A 28 2 2 2 2 2 2 2 2 2 2 2 2 2 3 2 2 2 ee fe he he he fe fe fe Ie 2c ie ie 2s Ie it 


void inputBinaryString(); 
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J [FR FF a 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 28 2 2 2 ee fe he fe fe fe he he ie ie Fie ie Ie ik 


// Method countRunningOnes 

// Return value int nclass, largest number of running ones 
//Parameters int start, int finish 

// Purpose counts the number of running ones in each necklace 
J [78 FR AR AR A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee fe he fe fe fe he fe Ie Ie ie ie iS Ie iS 


int countRunningOnes(int start, int finish); 


J [7 FR A A 2 2 2 2 2 2 2 2 He 2 2 2 2 2 2 2 2h 2 2 2 ee fe fe he he fe he he ie 2g Fie ie Ie ie 


// Method countRunningZeros 

// Return value int nclass, largest number of runnning zeros 

//Parameters int start, int finish 

// Purpose counts the largest number of running ones in each sequence 
J [78 FR A RA A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee fe he he fe he he fe Ie 2k 2g Fie ie Ie it 


int countRunningZeros(int start, int finish); 


J [FR A A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2h 2 2 2 ee fe fe fe fe he he he ie ie Fie ie Ie is 


// Method shiftSequenceToNecklace 

// Return value none 

//Parameters none 

// Purpose given an arbitrary n-tuple, it rotates the sequence until it generates 
necklace 

J [78 7 A FB A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee he he fe he fe fe fe ie ie ie ie 28 Ie iS 


void shiftSequenceToNecklace(); 


J [78 7 A A A A A 28 2s 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee fe he fe fe fe he fe ie ie ie 2k ie Ie is 

// Method generateShiftingOnesSequence 

// Return value none 

//Parameters long long neckdec, decimal equivalent of necklace 


// long long numclump, number of missings after necklace 

// long distance, distance between adjacent signposts 

// long necknum, necklace number associated with each signpost 
// Purpose generates shifting ones sequence to use as signposts 


// Note: Entries are stored in a matrix format that is read from left to right where the last 
// entry of the previous row precedes the first entry of the current row. This ensures I 

// have enough memory to collect enough signposts 

J [RF A A a 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2h 2 2 ee he he he fe fe fe he he ie ie ie ie Ie it 

void generateShiftingOnesSequence(long long necdec, long long clmpsz, long distance, 
long necknum); 


J [RR A A 2 2s 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2h 2 2 2 2 ee fe fe fe fe he he ie ie ie 28 2S it 


// Method getWeight 
// Return value void 
//Parameters int num 


71 


// Purpose obtain user weight input for shifting ones sigposts 
J [7B FR PA i 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee he fe he fe fe he he ie ie ie ie ie 2S it 


void getWeight(int num); 


J [7B FR A Ha 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2h 2 2 2 ee fe he fe he fe he he ie ie ie ie Ie is 


// Method cycleshiftSequence 

// Return value none 

//Parameters int index, index of largest sequence of running ones 

// Purpose cyclically rotate sequence until largest number of 

// running ones is at the leftmost position of the sequence 
J [78 FR A AB A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ae 2 ee fe fe fe fe fe he fe ie ie ie 2c 2s 2k is 


void cycleshiftSequence(int index); 


J [FR A a 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2h 2 2 2 2 he fe fe fe he he he ie Ie ie Fie ie Ie it 


// Method createInputString 

// Return value none 

//Parameters long long input 

// Purpose creates binary string version of user input decimal 

// to be matched to a subsection of the DeBruijn sequence 
J [7B FR AR AR AB A 2 2 2 2 2 2 He 2 2 2 2 2 2 2 2 2 3 2 ee he fe he he he he he Ie Ie ie ie ie Ie is 


void createInputString(long long input); 


J [78 FR AR AB A A 2 28 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 2 ae ee he fe he fe he he ie 2k ie ie 2s Ie it 


// Method createFirstNecklace 

// Return value none 

//Parameters long long input 

// Purpose creates first necklace to be used in creating a 
// subsection of the DeBruijn sequence 

J [7 FR A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2h 2 2 2 2 2 ee fe he fe he he he Ie ie ie ie ie Ie it 


void createFirstNecklace(long long input); 


J [78 A AB A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee he fe fe fe he fe ie 2k ie 2k ie Ie ie 


// Method createMiddleNecklace 

// Return value none 

//Parameters long long input 

// Purpose creates middle necklace to be used in creating a 
// subsection of the DeBruijn sequence 

J [7B FR A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2h 2 2h 2h 2 2 ee fe he fe fe he he he ie 26 ie ie Ie it 


void createMiddleNecklace(long long input); 


J [FR Ha A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee fe fe fe fe fe he Ie ie ie Fie ie Ie it 

// Method createLastNecklace 

// Return value none 

//Parameters long long input 

// Purpose creates last necklace to be used in creating a 
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// subsection of the DeBruijn sequence 
J [7 FR FR A Hi 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2h 2 2 2 2 he fe he fe fe fe he he ie ie ie ie Ie it 


void createLastNecklace(long long input); 


J [RFR FR A A 2 2 2 2 2 2 2 2 He 2 2 2 2 2 2h 2 2h 2 2 2 ee fe he fe fe fe he he ie ie ie ie Ie it 


// Method createDeBruijnSection 

// Return value none 

//Parameters none 

// Purpose concatenates nonperiodic portions of three adjacent necklaces 
J [7B FR A Hi a 2 2 2 2 2 2 2 2 2 2 2 2 2 2h 2 2h 2 2 2 ee fe he fe fe fe he he ie 2g ie 28 Ie is 


void createDeBruijnSection(); 


J [FR FR A Hs 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee fe he fe fe he he he ie ie ie ie Ie ie 


// Method createWraparoundDeBruijnSection 

// Return value none 

//Parameters none 

// Purpose concatenates nonperiodic portions of four adjacent necklaces: 
// LOO 0G FO 14 1110 


J [78 FR A AR A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 2 ee fe he fe he fe fe fe he 2k ie 26 ie 2k ie 


void createWraparoundDeBruijnSection(); 


J [78 FR FA 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2h 2 2h 2 2 ee he fe he fe fe he he he ie ie Fie ie Ie it 


// Method getInputDec 

// Return value none 

//Parameters long long input 

// Purpose obtains user input decimal value that yields binary string 
// to be searched for in DeBruijn sequence 

J [FR FR AR AR A A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee fe he he he he he he ie 2k ie ie it Ie is 


void getInputDec(long long input); 


J [7B AR AR A A 2 2 2 2g 2 2 2 2 2 2 2 2 2 2 2 2 3 2 ae 2 fe he he he fe fe he 2k 2k ie 2k ie Ie is 


// Method searchDeBruynString 

// Return value long long position 

//Parameters none 

// Purpose searches for user input binary string in DeBruijn sequence 
// and returns its location 

J [FR FR A Hs 2 2 2 2 2 2 2 2 He 2 2 2 2 2 2h 2 2h 2 2 2 2 2 fe he fe fe fe he he ie ie ie ie 2S is 


long long searchDeBruijnString(); 


J [78 FR AR A A A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee fe he fe he fe he fe he ie ie ie ie Ie is 


// Method findSequencePeriod 
// Return value int period 
//Parameters none 


// Purpose find period of a necklace 
J [FR FR A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2h 2 2 2 ee fe fe fe fe fe he he ie ie Fie ie Ie is 
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int findSequencePeriod(); 


J [7878 AR AB A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee he fe fe he fe fe fe ik 2k 2k ie ie Ie is 


// Method computeDensity 

// Return value int nclass, largest number of runnning zeros 

//Parameters int start, int finish 

// Purpose counts the total number of ones in the sequence 

// ** Note: To compute number of zeros = length of sequence - computeDensity 
J [FR A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee he fe he fe fe fe he he ie ie Fie ie IS it 


int computeDensity(int start, int finish); 


J [78 FR A AB A AB A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee fe he fe he fe he fe ie ie ie ie iS Ie is 


// Method testforNecklace 

// Return value bool foundnecklace 

//Parameters long long bininput 

// Purpose tests n-tuple to see if it is a necklace 
J [7 FR FR A i 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2h 2 2 2 ee he he fe he fe he he ie ie Fie ie Ie it 


bool testforNecklace(long long bininput); 


J [7B FR A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2h 2 2 ee he fe he fe he fe he Ie ie ie ie ie 2g ik 


// Method countNecklacesPerClass 

// Return value none 

//Parameters int ones, number of leading ones 

// Purpose separates the number of leading ones in each necklace 
// into distinct classes 

J [7 FR A A A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 2 ee he he fe fe fe he fe he 2k ie ie ie Ie is 


void countNecklacesPerClass(int ones); 


J [787 A A AB A 28 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee he he fe he he he fe he 2k ie ie 2s Ie is 
// Method generateChangeSequence 


// Return value none 
//Parameters long long oldec, decimal value of previous necklace 


// long long newdec, decimal value of current necklace 

// int indx, array index for sequence 

// Purpose counts the number of bits that differ between each necklace 
// and inserts these values into a sequence 


J [78 7 A A A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 2 ee he he fe fe fe he fe ie 2k 2g ie 2s Ie is 


void generateChangeSequence(long long olddec, long long newdec, int indx); 


J [78 78 A AB A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ae 2 he he he fe fe fe fe ie 2k ie ie ie Ie is 


// Method displayNecklaceInfo 
// Return value none 
//Parameters int col, classnumber, long long value, numclump 


// Purpose displays all information related to a necklace 
J [7B AR AR A A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 a 2 ee fe he fe he fe he fe ie 2k Ie 26 2s Ie is 
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void displayNecklaceInfo(int a, int b, long long c, long long d); 


J [78 7 A A A A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee fe he fe he fe he fe Ie 2k ie ie 28 Ie ik 


// Method displayPower2Array 

// Return value none 

//Parameters none 

// Purpose displays power2[i] array 


J [78 FR AR AR A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee fe he fe fe fe fe fe he 2k ie ie ie Ie is 


void displayPower2Array(); 


J [7B FR A AR A A 2 2 2 2 2 2 2 2 2 2 2 2g 2 2 2 2 3 2 ee he he fe fe fe he fe ee ie ie ie 2k ik 


// Method displayDifference 

// Return value none 

//Parameters none 

// Purpose calculates and displays difference between power2 approximation and 
number of necklaces 

J [7B FR A a 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2h 2 2 2 2 2 ee fe he fe he he he he ie ie Fie ie Ie it 


void displayDifference(); 


J [7B FR AR AR A A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee fe he fe he fe he he ie 2k IE ie iS IE is 


// Method displayNumbNeck 
// Return value none 
//Parameters none 


// Purpose displays number of necklaces per class 
J [FRR A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee fe he fe fe fe he ie ie ie Fie ie Ie it 


void displayNumbNeck(); 


J [78 FR ARR A A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 2 ee fe he he he he he fe ie Ie ie ie is Ie is 
// Method displayClassNumbers 

// Return value none 

//Parameters none 


// Purpose displays class numbers of necklaces in decreasing order 
J [FR RP A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee he fe fe fe he fe he ie ie ie ie ie 2S is 


void displayClassNumbers(); 


J [7B FR A AR A A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee fe fe he he fe he fe Ie Ie Ie ie ie Ie is 


// Method displayNcount 
// Return value none 
//Parameters none 


// Purpose displays total number of necklaces 
J [7 FR A A 2 2 2s 2 2 2 2 2 2 He 2 2 2 2 2 2 2 2 2 ee he he he fe fe fe he he ie ie Fie ie Ie it 


void displayNcount(); 


J [7B FR A AR A A A 2 2 2 2 2 He 2 2 2 2 2 2 2 2 2 3 2 ee he he fe he he he he Ie Ie iC ie ie IE iS 


// Method displayClumpsPerClass 
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// Return value none 

//Parameters none 

// Purpose displays total number of missing even n-tuples per class 
J [7B A AR A AB A 2 2 2 2 2 He 2 2 2 2 2 2 2 2 2 3 2 ee he he he he fe he fe Ie 2k ie 2k ie Ie ik 


void displayClumpsPerClass(); 


[[PEERERERE ERB BER RRR RE EG ee eR a RG Re 


// Method displayClumpDistribution 

// Return value none 

//Parameters none 

// Purpose displays various clump sizes per class 
J [7 FR A AR A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 2 2 ee he fe he fe he he ie Ie Ie ie ie Ie iS 


void displayClumpDistribution(); 


J [7 FR A AR A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 2 ee he he he he he he he Ie Ik ie ie iC IE is 


// Method displayDecimalArray 

// Return value none 

//Parameters none 

// Purpose displays decimal equivalent of necklaces associated with clumps per 
class 

J [FR FR A A 2 2s 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2h 2 2 2 ee fe he fe fe he he Ie ie ie ie ie 2S it 


void displayDecimalArray(); 


J [78 FR A AB A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee he fe fe he fe he fe 2k Ie ie 26 2s Ie is 


// Method displayShiftingOnesSignposts 

// Return value none 

//Parameters none 

// Purpose displays decimal and binary equivalents of shifting signposts 
J [78 FR A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee fe he fe fe he he Ie ie ie Fie ie Ie it 


void displayShiftingOnesSignposts(); 


J [7 A AS A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee he he fe he fe he fe ik 2k ie ie ie Ie is 


// Method displayShiftingOnesDistances 

// Return value none 

//Parameters none 

// Purpose displays distances between shifting ones signposts 
J [78 FR A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2h 2 2 2 ee fe he fe fe fe he he ie ie Fie ie Ie it 


void displayShiftingOnesDistances(); 


J [78 A AR A A 2 2 2 2 2 2 He 2 2 2 2 2 2 2 2 2 2 2 ee fe fe fe he fe he fe he 2k 2k ie Ie Ie iS 


// Method display BitChangeSequence 
// Return value none 
//Parameters none 


// Purpose displays bit changes between necklaces 
J [7 RR A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee fe he fe fe fe he Ie ie ie Fie ie 2S is 
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void displayBitChangeSequence(); 


J [78 FR A AS A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 2 ae 2 fe he fe he fe he fe he 2k Ie 2c ie Ie ik 


// Method displayFirstNecklace 

// Return value none 

//Parameters none 

// Purpose displays prior necklace in deBruijn section 
J [78 AR AB A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 2 ee fe he fe he fe fe fe ie 2k ie ie ie Ie ie 


void displayFirstNecklace(); 


J [7 AR AB A A 2 2 2s 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee fe he fe he he he fe ie 2k ie ie ie Ie is 


// Method displayMiddleNecklace 
// Return value none 
//Parameters none 


// Purpose displays input necklace in deBruijn section 
J [78 FR A AR A A 2 A 2 2g 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee he he fe fe fe he fe he 2k ie ie ie Ie is 


void displayMiddleNecklace(); 


J [7B 8 A AB A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 2 he he fe he fe he fe ie ie ie ie ie ie is 
// Method displayLastNecklace 

// Return value none 

//Parameters none 


// Purpose displays post necklace in deBruijn section 
J [77 AR AB A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 2 ee fe he fe he he fe fe he ie ie 26 ie Ie is 


void displayLastNecklace(); 


J [7 FR A AR A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ae ee he fe fe fe he fe he 2k 2c ie ie Ie is 
// Method displaydeBruijnSection 

// Return value none 

//Parameters none 


// Purpose displays deBruijn sectional associated with user input necklace 
J [78 7 A AR A A 2 2 2 2 2 2 2 2 He 2 2 2 2 2 2 2 2 2 ae 2 he he fe fe fe fe fe 2k 2k ie ie 28 Ie is 


void displaydeBruijnSection(); 


J [7B A AB A A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 2 3 2 he he fe fe fe he fe ie 2k 2k ie ie Ie is 


// Method display UnobservedSequences 
// Return value none 
//Parameters none 


// Purpose displays all missings (even binary n-tuples) between any 2 nunmbers 
J [7B 7 AR AR A A 2 2 2 2 2 2 He 2 2 2 2 2 2 2 2 2 2 2 ee ee fe he fe he fe he 2k 2g ie ie Ie is 


void display UnobservedSequences(); 


J [78 AR AB A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee he he fe fe fe fe he ie ie ie 2k 2 2k 


// Method outputData 
// Return value none 
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// Parameters 
// Purpose display numerical results of algorithm output 
J [78 FR FR AB A A 2 2 2 2 2 2 2 2 2 2 2 2 2 3 2 3 ee fe he he fe fe he he ie ie ie 2k ie 2k 


void outputDataQ); 


private: 


J [FR FF A 2 2 2s 2 2 2 2 2 2 2 2 2 2 2 2 2 2h 2 2 ee he fe he fe fe fe he he ie ie Fie ie Ie it 


// Method initSequences 

// Return value none 

//Parameters int size, int thresh, float per 

// Purpose to initialize all the sequences involved in the analysis of the algorithm 
J [7B FR A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2h 2 2 2 ee fe he fe he he he he ie ie Fie ie Ie it 


void initSequences(int size, int thresh, int per); 
private: 


//variables 

int m; // size of sequence 

int k; // number of leading ones input by the user 

int q; //index for onescount array 

int mark; // marker for shifting ones array 

int dim; //dimension for decimalarray 

int begin; //index to begin counting ones or zeros 

int threshold; // smallest size of missings displayed 

int num; // integral numerator of percentage 

int numevensignposts; // number of evenly spaced signposts 

int ncountdiff; // counts the number of necklaces between adjacent signposts 

int firstperiod; // period of first adjacent necklace in a deBruijn section 

int middleperiod; // period of middle adjacent necklace in a deBruijn section 

int lastperiod; // period of last adjacent necklace in a deBruijn section 

int sumperiod; // sum of last three periods 

int weight; // weight for shifting ones signposts 

int rowindex; // shifting ones row index <= MAX_LENGTH 

int colindex; // shifting ones column index <= MAX_COLUMNS 

long distance; // computes distance between initial and final necklace 

long olddistance; // computes distance of last necklace 

long totaldistance; //distance from end of deBruijn sequence (111.11) until first 
//necklace in deBruijn section 

float percent; // percentage of maximum clump size in a given class 

long ncount; // necklace number on theta generated list 

long noncount; // theta-generated non-necklace count 

long shiftonesnoclump; //number of shifting ones sequence NOT associated with 

//missings 
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long long absposition; // // absolute position from origin or end of 000...0 necklace 
long long totmissing; // number of missing even n-tuples in a given class 

long long maxclumpsize; // maximum clump size for all the classes 

long long inputstrdec; // decimal equivalent of input string 

double shiftonesclump; //number of shifting ones sequence associated with missings 
double binlength; // binary length of 2m-1 

double petriumilestones; // number of petriu milsetones assuming max period of m 
double milestonedistance; // petriu distance given number of signposts 

double topt; // optimal petriudistance between milestones 


//arrays 
int sequence[MAX_LENGTH]; // binary sequence for necklaces 
int divisorssMAX_LENGTH+1]; // all divisors for a given n 
int bitchangeseqiMAX_COLUMNS]; // stores the number of bit changes between 
// necklaces 
int inputstring;j MAX_LENGTH]; // binary sequence of user input decimal 
int firstnecklace[~MAX_LENGTH]; // necklace listed before user input necklace 
int middlenecklace[~MAX_LENGTH]; // user input necklace 
int lastnecklace[MAX_LENGTH]; // necklace listed after user input necklace 
int debruijnsectionIMAX_LENGTH]; // section of DB sequence containing prior 
// three arrays 
int wrapdebruijnsection[MAX_LENGTH]; // section of DB sequence that includes 
// strings between 00...0 and 11...1 
long onescountt{MAX_LENGTH+1]; //collects number of necklaces per class 
long power2;| MAX_LENGTH+1]; // power of 2 approximation to onescount 
long diff[MAX_LENGTH+1]; // difference between power2 and onescount 
long ncountarray[MAX_LENGTH][MAX_COLUMNS]; // numbered count of each 
//necklace on theta list per class 
/Aong long decimalarray[MAX_LENGTH][MAX_COLUMNS]; // decimal 
//equivalent of each necklace per class 
long clumpcntarray[MAX_LENGTH+1]; //number of clumps per class 
long long missingarray[MAX_LENGTH][MAX_COLUMNS]; 
long long maxclumpsizearray[MAX_LENGTH]; // maximum clump size for each 
// class 
long long totalmissings[|MAX_LENGTH]; 
long long missingbelowthreshold[MAX_LENGTH]; // total number of missing for a 
// given class less than threshold 
long long clumpsizearray[MAX_LENGTH][MAX_ COLUMNS]; 
long long evenspacedsignpostarray[MAX_LENGTH][MAX_ COLUMNS]; //contains 
// decimal value of evenly spaced signposts 
long long shiftonesarray[MAX_LENGTH][MAX COLUMNS]; // contains decimal 
//value of necklaces with shifting ones in a given class 
long shiftonesdistancearray[MAX_LENGTH][MAX_COLUMNS]; // contains 
// distances (in bits) between adjacent signposts 
long shiftonesncountarray[MAX_LENGTH][MAX_COLUMNS]; // contains 


79 


// necklace number of each signpost 
long long shiftonesclumpsizearray[MAX_LENGTH][MAX_COLUMNS]; // contains 
// clumpsize associated with shifting ones array 
}; //end Class Necklaces 


#endif // end of file Necklaces.h 


2. IO Thesis, FileOpeningException and SieveSizeException Header 
Files 


J [5 7 2 2 2g 2g ae 2 ae ae a 2 ee ee fe he he he he he ie 26 28 26 28 28 28 28 2s 28 2 2s 2s 22s 2 2 2 2h 2 2h 2h 2h 2h 2 2 2 2 2 2 2 2 2 2 he fe fe fe he fe he he 2c ie ie ie 2s it 
oh K ok 


// 1OThesis.h 

// 

// LT John Ortiz 

// 

// 

// Course: CS2971 2006 QI 

// Project: Generating Necklaces 
// Operating Environment: Windows XP Home 
// Compiler: Visual Studio .NET 
// Date: 08 December 2005 

// Description: 

// 


J [7 7 2 2 2 2g ae ae ae ae ae ae he he he he he he fe he he 6 26 26 28 28 28 28 28 2h 28 28 28 28 2s 22 2 2 2 2 2 2k 2 2h 2 2 2 28 2 2 2 3 2 ee fe fe fe he fe fe fe fe Ie ie ie ie Ie ik 
ok Kk 


#ifndef IOTHESIS_H 
#define IOTHESIS_H 


#include <string> 

#include <sstream> 
#include "NecklaceTable.h" 
#include <fstream> 


using std::string; 
using std::ifstream; 


const int EXIT_WITH_ERROR = 1; // global constant that enables a program to exit 
when there is an error. 

const int EXIT_NORMALLY = 0; // global constant that enables a program to exit when 
there is no error. 


[7 7 2 2 2g ae ae ae ae ae ae 2 ae he fe fe fe fe he he ie ie ie 26 28 ie 2s 28 2s 2s 2s 2s 2s 2 2 2 2 2 2 2k 


// Class: JTOThesis 
// Purpose: Display a 20x20 Matrix of integers where multiples of a filter 
// provided by the user are replaced by blank spaces. 
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class IOThesis { 


public: 
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// Method displayMessage 

// Return value none 

//Parameters none 

// Purpose prompts user for dimension input 
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static void displayMessage(string message); 
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// Method openOutputFile() 

// Return value 

// Parameters: const char *FILE_NAME the name of the file to open for writing 
// Purpose: dynamically create a new ofstream and returns a pointer if successful 
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static ofstream* openOutputFile(const char *FILE_NAME); 
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// Method openInputFile 

// Return value ifstream* 

// Parameters const char *FILE_NAME, the file to open for reading 
// Purpose 
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static ifstream* openInputFile(const char *FILE_NAME); 
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// Method pauseBeforeExit 

// Return value none 

//Parameters int condition (EXIT_WITH_ERROR = 1, global defined in 
IOController.h) 

// (EXIT_NORMALLY = 0, global defined in IOController.h) 

// Purpose keeps display window open long enough to see results 
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static void pauseBeforeExit(int condition); 
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// Method inputInt 

// Return value int 

//Parameters none, function will ask for input via I/O 

// Purpose take string input and determine if valid for an int 


// use inputFilter = [OController::inputInt(); vice cin >> inputFilter 
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static int inputInt(); 
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// Method inputIntWithLimits 

// Return value int 

//Parameters (string message, int lower, int upper) 

// Purpose calls input int for a valid int, then checks limits 
// used similiarly to inputInt 
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static int inputIntWithLimits(string message, int lower, int upper); 
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// Method clearInputBuffer 

// Return value bool TRUE if cin buffer has no additional stuff in it 
//Parameters none 

// Purpose utility function for processing inputs 
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static bool clearInputBufferQ; 
}; //end Class IOThesis 


#endif // end of file I[OThesis.h 
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// FileOpeningException.h 

// 

// LT John Ortiz 

// 

// 

// Course: CS2971 2006 QI 

// Project: 

// Operating Environment: Windows XP Pro 
// Compiler: Visual Studio .NET 

// Date: 08 December 2005 


// Description: Provide exception handling for file opening. 
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#include <stdexcept> 
using std::runtime_error; 
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// Class: FileOpeningException 


// Purpose: Provide exception handling for file opening. 
[ [PRR RRR A EE RR ER EE eae EE EE RR EE RE 


class FileOpeningException : public runtime_error { 
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public: 
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// Method FileOpeningException--Constructor 
// Return value none 
//Parameters none 
// Purpose 
J [78 FR A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 a 2 he fe he fe fe fe he ie ie ie ie ie 2s 
FileOpeningException::FileOpeningException(): runtime_error("Unable to open file") 
{ 
//nothing else to initialize 


} 


//variables 
private: 
//none 
}; //end Class FileOpeningException 
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// SieveSizeException.h 

// 

// LT Ortiz 

// 

// 

// Project: Thesis 

// Operating Environment: Windows XP Home 
// Compiler: Visual Studio .NET 


// Date: 

// Description: Provide exception handling for incorrect necklace size requests 
// There is no default constructor. Coders must provide an argument. 
// Choose the appropriate enumerated constant 
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#ifndef SIEVESIZEEXCEPTION_H 
#define SIEVESIZEEXCEPTION_H 


#include <string> 
#include <stdexcept> 
using namespace std; 
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// Class: SieveSizeException 


// Purpose: Provide exception handling for incorrect sieve size requests 
J [7878 AR AR A A 2 2 2 2s 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee he he fe he fe he fe ie Ie Ie ie 28 2k it 


83 


class SieveSizeException : public runtime_error { 


public: 
const static enum { TOO_SMALL, TOO_BIG }; 
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// Method SieveSizeException--Constructor 

// Return value none 

//Parameters const int ERROR_ID 

// Purpose 
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SieveSizeException(const int ERROR_ID) : runtime_error( 
resolveErrorType(ERROR_ID) ) { 

//nothing else to initialize 


} 
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// Method resolveErrorType 

// Return value string The exception message 
//Parameters const int errorID 

// Purpose generate exception text 
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string resolveErrorType(const int ERROR_ID) { 


switch (ERROR_ID) { 
case TOO_SMALL: 
return "A necklace needs to have positive dimensionality"; 
break; 


case TOO_BIG: 
return "That would produce a necklace too large to output in the console"; 
break; 


default: 
return "Really bad things are happening here...function was not provided 
TOO_SMALL or TOO_BIG as an argument"; 


} 


}//end function resolveErrorType 


//variables 
private: 
//none 
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}; end Class SieveSizeException 


#endif // end of file SieveSizeException.h 
oF CPP FILES 
L. Necklace CPP File 
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// Necklaces.cpp 

// 

// LT John Ortiz 

// 

// 

// Project: Thesis 

// Operating Environment: Windows XP Home 
// Compiler: Visual Studio .NET 
// Date: 

// Description: 

// 


[28987 RFR RE A A EE 2 A EA 2 2 C2 2A EA EO 2 2 2 E22 2 ICH HAE 2 2 CH EH oO OE 
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#include <iostream> 
#include <iomanip> 
#include <cctype> 
#include "math.h" 
#include "Necklaces.h" 
#include "IOThesis.h" 


using namespace std; 
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// Method Necklaces--Default Constructor 
// Return value none 
// Parameters none 
// Purpose 
J [7B A AB A A 2 2s 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee he fe he he fe he he ie Ie ie 2k ie is 
Necklaces::Necklaces() 
{ 
initSequences()MAX_COLUMNS,DEFAULT_THRESHOLD,DEFAULT_PERCENT); 
/1/ MAX COLUMNS = 34 
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// Method Necklaces-- Constructor 
// Return value none 
// Parameters int size, int thresh, float per 
// Purpose 
J [7B FR FR A a 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee fe he fe fe fe he he ie Fie ie ie Is 
Necklaces::Necklaces(int size, int thresh, int per) 
{ 

initSequences(size, thresh, per); 
}//end constructor 
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// Method initSequences 

// Return value none 

//Parameters int size, int thresh, float per 

// Purpose to initialize all the sequences involved in the analysis of the algorithm 
[7 7 2 2 ag 2g 2g ae ae ae ae ae ae he he fe fe fe he he ie ie he 26 28 28 2s 2s 2s 2s 2s 2s 2s 2 2 2 2 2 2 2k 


void Necklaces::initSequences(int size, int thresh, int per) 


{ 
m = size; // defined only once for all functions 
threshold = thresh; 
begin = 0; 
num = per; 
percent = num /100; 


if(m<1) 
{ 
cout << ''A necklace needs to have a length > = 1 \n" 
<< "exiting the program.\n"' << endl; 


IOThesis::pauseBeforeExit(EXIT_WITH_ERROR); 
} 


if (m > MAX_LENGTH) // MAX_COLUMNS = 400 
{ 


cout << ''That would produce a necklace too large to output in the console \n" 
<< "exiting the program.\n" << endl; 


IOThesis::pauseBeforeExit(EXIT_WITH_ERROR); 
} 


86 


initNecklace(0); //initializes necklace to be all zeros 


for (int i = 0; 1 << MAX_LENGTH; i++) 
{ 
maxclumpsizearray[1i]= 0; 
totalmissings[i] = 0; 
missingbelowthreshold[i] = 0; 
firstnecklace[i] = 9; 
middlenecklace[i] = 9; 
lastnecklace[i] = 9; 
debruijnsection[i]= 9; 
wrapdebruijnsection[i] = 9; 


} 
for (int i = 0; 1 <= MAX_LENGTH,; i++) 
{ 

power2[i] = 1; 

diffi] = 0; 


onescount[i] = 1; 
clumpentarray[i] = 0; 
divisors[i] = -1; 


} 


for (int i= 0; 1 << MAX_COLUMNS; i++) 
{ 

bitchangesegq[i] = 0; 
} 


mark = -1; 
for (int i= 0;1< MAX _ LENGTH; i++) 
{ 
for(int j = 0; } < MAX_COLUMNS; j++) 
{ 
missingarray[1][j] = 0; 
ncountarray[i][j] = 0; 
//decimalarray [i][j]= 9; // a necklace will always have an even decimal equivalent 
clumpsizearray[1][j] = 0; 
shiftonesarray[1i][j] = mark; 
shiftonesdistancearray[i][j] = 0; 
shiftonesclumpsizearray[1i][j] = 0; 
evenspacedsignpostarray[i][j] = mark; 
} 
} 


for (int i= 0; 1 < MAX LENGTH; i++) 
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{ 
for(int j = 0; j < MAX_LENGTH; j++) 


evenspacedsignpostarray[i][j] = 0; 
} 
} 


} // end initSequences 
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// Method initNecklace 

// Return value decimal equivalent 

// Parameters int ones 

// Purpose to create a sequence consisting of a certain number of leading ones 
// with the remaining sequence consisting of zeros 
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long long Necklaces::initNecklace(int ones) 


{ 


long long initdec = 0; 


for (int i=0; 1 < m; i++) 
{ 


sequence[i] = 0; // sequence initially set to all zeros 


} 


k = ones; 
int x = 1; 
for(int i=0; 1< k; i++) // initializing initial fill with given number of leading ones 
{ 
sequence[i]= x; 
} 
initdec = binToDec(); 


return initdec; //returns decimal equivalent of initial necklace 
}// end initNecklace 
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// Method printSequence 

// Return value none 

// Parameters none 


// Purpose displays binary sequences into three parts 
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void Necklaces::printSequence() 


{ 


int n = 0; 
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inti=0; 
while(sequence[i] != 0 && i< _m) // determine number of leading ones 


{ 
n++; 
I++; 
} 
for(int i = 0; 1 < m; i++) 
{ 
if(i == n) 
{ 


wom, 
bf 


cout << sequence[i] << 

// this next line is optional and can be commented out if desired 

//cout << setw(OUTPUT_CELL_SIZE - 1)<<""; //print 2 spaces after first run of 
ones and a zero 


else if(i == m-2) 

{ 
cout << sequence[i] << 
// this next line is optional and can be commented out if desired 
//cout << setw(OUTPUT_CELL_SIZE - 1)<<"'"; //print 2 spaces before 

terminating zero 

} 

else 

{ 
cout << sequence[i] <<" "; //print 1 within middle part 

} 

} 


}// end printSequence 


mom, 
2 
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// Method runNecklaceAlgorithm 
// Return value none 
//Parameters long long initial, long long final, long long input, bool numflag 
// Purpose executes necklace algorithm 
J [RFR FR A A 2 2 2 2 2 2 2 2 ae He 2 2 2 2 2 2 2h 2 2 2 2 2c fe fe fe he fe he he ie ie ie ie Ie is 
void Necklaces::runNecklaceAlgorithm(long long initial, long long final, long long input, 
bool numflag) 
//void Necklaces::runNecklaceAlgorithm(int ones, long long numneck, bool numflag) 
{ 
ncount = 0; 
noncount = 0); 
totmissing = 0); 
maxclumpsize = 0; 
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distance = 0); 

olddistance = 0; 

ncountdiff = 0; 

shiftonesclump = 0; 

shiftonesnoclump = 0; 

firstperiod = 0; 

middleperiod = 0; 

lastperiod = 0; 

sumperiod = 0; 

int classnumber = 0; 

int j = m-1; //largest index such that aj=1 and ak=0 for k=j+1 to m 

//k = ones; 

q=0; 

long inputncount = 0; // necklace number associated with user input necklace 

long long nondecimal = 0; // decimal equivalent of current non-necklacel 

long long ndecimal = 0; // decimal equivalent of current necklace 

long long oldndecimal = 0; // decimal equivalent of previous necklace 

long long oldnondecimal = 0; // decimal equivalent of previous non-necklace 

long long priorndecimal = 0; 

long long maxnum = final - initial + 1; 

long long value = 0; 

long long inputdec = input; 

long long nummissing = 0; //difference between decimal equivalences of 2 adjacent 
necklaces 

int col = 0; // column number for arrays 

int idx =0; // index number 

int idx2 = 0; // second index number 

int oldclassnumber = 0; 

int numclump = 0; // number of clumps in a given class corresponding to number of 
signposts 

int numevensignposts = 0; // number of evenly spaced signposts 

bool printflag = numflag; 

rowindex = 0; 

colindex = 0; 


//initNecklace(k); 


ncount = |; 

shiftonesnoclump++; 

oldndecimal = initial; 

decToBin(oldndecimal); // initial necklace 

distance = findSequencePeriod(); // initial calculation 


if (printflag == true) // print statement for initial necklace 


{ 
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cout << endl; 

cout << ''\nNecklace '' << ncount << '":'' << endl; 
printSequence(); 

cout << '' decimal equivalent = '' << oldndecimal << end]; 
cout << "\n"; 


} 


do 
{ 
nummissing = 0; //resetting the count 
j=m-1; //always start at the end of the current sequence and work backwards to find 
largest aj = 1 
while (sequence[j]==0) // keep decrementing index until aj = 1 


jai-1; 
} 


sequence[j]=sequence[j]-1; //subtracts | from jth position 


ifG!=m-1)// checks to see if we can copy beginning portion to end 
{ 
for(int i=1; 1 < m -j; i++) 
{ 
sequence[j+i] = sequence[i-1]; //copies the first portion onto the remainder of 
necklace 


} 
} // end if 


if(m %(j+1)!=0) // Checking for generated nonnecklaces: (j+1) ensures we do not 
have division by zero; m = real length, j+1 = real index 
{ 
noncount = noncount++; 
nondecimal = binToDec(); 
/Af (printflag == true) 


//cout << endl; 

//cout << "\nNon-Necklace " << noncount << ":" << endl; 
//printSequence(); 

//cout << "decimal equivalent = " << nondecimal << endl; 
//cout << "\n" ; 


} 
} // end if 


if(m %(j+1)==0) // Checking for necklaces: (j+1) ensures we do not have division by 
zero 


{ 
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ncount = ncount++; 

ndecimal = binToDec(); 

distance += findSequencePeriod(); // computing and summing up successive 
// periods 

olddistance = distance - findSequencePeriod(); 


if(count == maxnum) // allows us to exit after a given number of thetasteps 
{ 
break; 


} 


if(mdecimal == inputdec) // case where missing n-tuple is close to its necklace 
{ 
//cout << "\noldndecimal for first necklace = " << oldndecimal << endl; 
createFirstNecklace(oldndecimal); 
firstperiod = findSequencePeriod(); 
//cout << "\nndecimal for middle necklace = " << ndecimal << endl; 
createMiddleNecklace(ndecimal); 
middleperiod = findSequencePeriod(); 
inputncount = ncount; 
totaldistance = olddistance - firstperiod; // need to adjust for first necklace 


} 


if(final < inputdec) 


{ 
if(ncount == inputncount + 1) 
{ 
//cout << "\nndecimal for last necklace =" << ndecimal << endl; 
createLastNecklace(ndecimal); 
lastperiod = findSequencePeriod(); 
sumperiod = firstperiod + middleperiod + lastperiod; 
} 
} 


if Gdx <= MAX_COLUMNS) // prevents access violations in memory 
{ 
generateChangeSequence(value, ndecimal, idx); 
idx++; 
} 
if(computeDensity(0, m-1) <= weight + classnumber)// tests for shifting ones in 
substring: changed classnumber to 0 
{ 
//generateShiftingOnesSequence(ndecimal, nummissing, distance, ncount); 
shiftonesnoclump++; 
} // end if(computeDensity) 
//{ 
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// evenspacedsignpostarray[col][col] = ndecimal; 

//} 

if (printflag == true) 

{ 
cout << endl; 
cout << '\nNecklace "' << ncount << '':" << endl; 
cout << '\n"; 
printSequence(); 
cout << ''"" decimal equivalent = '"' << ndecimal << endl; 
//cout << "\nmaxdistance = " << distance << endl; 
cout << '\n"'; 

} // end if(printflag) 

classnumber = countRunningOnes(begin,j); 

countNecklacesPerClass(classnumber); 

nummissing = (oldndecimal - ndecimal - 2)/2; 

if ((nummissing) > 0) // indicates there exists a clump 


value = 0; //resetting value to accept current value 
value = oldndecimal; // sets "value" to previous value of ndecimal 
if (printflag == true) 
{ 
cout << "number missing between Necklaces '' << ncount - | <<" and "" << 
ncount << " =" << nummissing << endl; 
cout << '"\n"; 
} // end if(printflag) 
totmissing += nummissing; 


if(oldclassnumber != classnumber) // verifies if class number has changed 


oldclassnumber = classnumber; 
col = 0; 
numclump = 0; // reset when class number changes 
/fidx = 0; // check to see if you need it here 
} 
else 
{ 
col++; // change to allow new info in a given class to be recorded 
totalmissings[classnumber]+= nummissing; // sums total missing for a given 
class 
} // end if-else 
if (Nummissing > maxclumpsize) // only records largest clump size per class 
{ 
maxclumpsize = nummissing; 
maxclumpsizearray[classnumber]= maxclumpsize; 


} 
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if(nummissing >= threshold && col <= MAX_COLUMNS) // restricting to 
collecting signpost information 
{ 
clumpsizearray[classnumber][col] = nummissing; 
//decimalarray[classnumber][col] = value; //stores decimal equivalent of 
sequence after missings 
//cout << "decimalarray[" << classnumber << "][" << col <<"] =" << value; 
//cout << "\n" ; 
missingarray[classnumber][col] = nummissing; // records instances of all 
missings > threshold for a given class 
idx++; 
} 
if(nummissing <= threshold) 
{ 
missingbelowthreshold[classnumber]+= nummissing; // records only the sum of 
missings below threshold 
} 
if(nummissing >= threshold) // counts the total number of clumps for a given class 
{ 
numclump++; 
clumpcentarray[classnumber]= numclump; 


} 


if(computeDensity(0, m-1) <= weight + classnumber)// tests for shifting ones in 
substring: changed classnumber to 0 


{ 


generateShiftingOnesSequence(ndecimal, nummissing, olddistance, ncount); 
shiftonesclump++; 

} // end if(computeDensity) 

ncountarray[classnumber][col] = ncount - 1; // why do I have this here?? 

if (printflag == true) 


//displayNecklaceInfo(col, classnumber, value, numclump); 


}// end outside if(nummissing) 
oldndecimal = ndecimal; // setting value which is the old "value" of ndecimal to 
current ndecimal 
///value = 0; //resetting value to accept current value 
} // end if(m %(j+1)==0) 
}while (ndecimal != final); // end of necklace algorithm 


}// end runNecklaceAlgorithm 


J [78 FR A AR A A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 2 ae 2 he he fe fe fe he he ie Ie Ie 28 Zit Ie is 
// Method inputBinaryString 
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// Return value none 
// Parameters none 
// Purpose checks whether user input binary string is valid 
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void Necklaces: :inputBinaryString() 
{ 
const int SIZE = MAX_LENGTH; 
char ch[SIZE]; // temporary buffer for input characters 
for(int i=0; 1i< SIZE; i++) //initializing ch array 
{ 
ch[i] = 9; 
} 
bool state = true; 
int limit = m; // limits size of user input binary string to match with constructor 


while(state == true) // This loop prompts a user to enter in only one integer that is less 
than 12 characters long 
{ 
intq =0; 
int total = 0; 
char element = cin.get(); 
while(element != '\n')// reads elements into a character array 


{ 
ch[q] = element; 
total++; 
if (total == limit) // checks if total exceeds 12 characters 
{ 


cout << 'You have reached the maximum number of characters allowed" << 
endl; 
cout << "We will truncate the additional digits" << endl; 
state = false; 
break; 
} 


if (cin.eof()) // checks if ctrl-z in input as a character 


cout << '\nInvalid character."; 
cin.clear(); //resets input stream so it can keep processing 


cin.putback('\n'); 

state = !TOThesis::clearInputBuffer(); 

break; 
} 
if (ch[q] != '0' Il ch[q] !='1') // checks if there is an invalid element 
{ 


cout << '\nInvalid character."; 
state = !1OThesis::clearInputBuffer(); 
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break; 


} 
else 
{ 
element = cin.get(); // gets the next element 
qtt; 
} 
} // end inside while 


} // end outside while 


for (int i=0; i <= m-1; i++) 


{ 
if(chfi]!= 9) 
{ 
sequence[i] = ch[i]; 
} 
} 


} // end inputBinaryString 
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// Method countRunningOnes 

// Return value int numones, largest number of runnning ones 
// Parameters int start, int finish 


// Purpose counts the largest number of running ones in each sequence 
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int Necklaces::countRunningOnes(int start, int finish) 
{ 

int alpha = start; 

int omega = finish; 

int tempent = 0; 

int numones = 0; 


for(int 1 = alpha; 1 <= omega; 1++) 


{ 
if(sequence[i] != 0) 
{ 
tempcnt++; 
} 
else 
{ 
tempcnt = 0; //resets tempcnt 
} 
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if(tempcnt > numones) 


{ 


numones = tempcnt; //ensures largest string of ones is preserved 


} 


} 
if(alpha == omega) 
{ 
numones = sequence[alpha]; 
} 


return numones; 
} // end countRunningOnes 
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// Method countRunningZeros 
// Return value int numzeros, largest number of runnning zeros 
// Parameters int start, int finish 


// Purpose counts the largest number of running ones in each sequence 
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int Necklaces::countRunningZeros(int start, int finish) 


{ 
int alpha = start; 
int omega = finish; 
int tempent = 0; 
int numzeros = 0; 


for(int 1 = alpha; 1 <= omega; 1++) 
{ 


if(sequence[i] != 1) 


tempcnt++; 


} 


else 


{ 


tempcnt = 0; //resets tempcnt 


} 


if(tempcnt > numzeros) 


{ 


numzeros = tempcnt; //ensures largest string of zeros is preserved 


} 
} 
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if(alpha == omega) 


{ 
if(sequence[alpha]== '0") 
{ 
numzeros = |; 
} 
else 
{ 
numzeros = 0; 
} 
} 


return numzeros; 
}// end countRunningZeros 
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// Method cycleshiftSequence 
// Return value none 
//Parameters int index, index of largest sequence of running ones 


// Purpose cyclically rotate sequence a number of positions to the left 
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void Necklaces::cycleshiftSequence(int index) 
{ 
int temp = 0; // temporary place holder 
int seqindx = index; 
for (int i = 0; i < seqindx; i++) //outer for loop only keeps track of number of shifts 
{ 
temp = sequence[0]; // beginning of sequence stored until it can be placed at the end 
for (int j = 0; j < m-1; j++) 
{ 
sequence[j] = sequence[j+1]; // all bits shift one place to the left 


} 


sequence[m- | ]= temp; 


} // end cycleshiftSequence 
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// Method createInputString 

// Return value none 

//Parameters long long input 

// Purpose creates binary string version of user input decimal 
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// to be matched to a subsection of the DeBruijn sequence 
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void Necklaces::createInputString(long long input) 


{ 


long long inputdecimal = input; 
decToBin(inputdecimal); 


for(int 1 = 0; 1 < m; i++) 
{ 
inputstring[i] = sequence[i]; 


} 


} // end createInputString 
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// Method createFirstNecklace 

// Return value none 

//Parameters long long input 

// Purpose creates first necklace to be used in creating a 
// subsection of the DeBruijn sequence 
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void Necklaces::createFirstNecklace(long long input) 
long long firstdecimal = input; 
decToBin(firstdecimal); 


for(int i = 0; 1 < m; i++) 
{ 
firstnecklace[i] = sequence[1]; 


} 


} // end createFirstNecklace 
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// Method createMiddleNecklace 

// Return value none 

//Parameters long long input 

// Purpose creates middle necklace to be used in creating a 
// subsection of the DeBruijn sequence 
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void Necklaces::createMiddleNecklace(long long input) 

t 


long long middledecimal = input; 
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decToBin(middledecimal); 


for(int i = 0; 1 < m; i++) 
{ 
middlenecklace[i] = sequence[i]; 
} 
} // end createMiddleNecklace 
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// Method createLastNecklace 

// Return value none 

//Parameters long long input 

// Purpose creates last necklace to be used in creating a 
// subsection of the DeBruijn sequence 
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void Necklaces::createLastNecklace(long long input) 
long long lastdecimal = input; 
decToBin(lastdecimal); 


for(int i = 0; 1 < m; i++) 
{ 
lastnecklace[i] = sequence[i]; 


} 


} // end createLastNecklace 
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// Method createDeBruijnSection 
// Return value none 
//Parameters none 
// Purpose concatenates nonperiodic portions of three adjacent necklaces 
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void Necklaces::createDeBruijnSection() 
{ 
int pl = firstperiod; // period of firstnecklace 
int p2 = middleperiod; // period of middlenecklace 
int p3 = lastperiod; // period of lastnecklace 
int p4 = 0; // period of last necklace during wraparound case 
long long testinteger = 0; //testing for wraparound case 


for(int 1 = 0; 1< pl; i++) 


{ 
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debruijnsection[i] = firstnecklace[i]; // first third of deBruijn section 
//cout << "\ndebruijnsection[" << i << "] =" << debruijnsection[i] << endl; 
//cout << "\nfirstnecklace[" << 1<<"] =" << firstnecklace[i] << endl; 


} 


//testinteger = binToDec(); 

//cout << "\ntestinteger = " << testinteger << endl; 

/Af(testinteger == generatePower2(m-1) Il testinteger == 0) // wraparound case 
INK 

// initNecklace(0); // all zeros 

// middleperiod = findSequencePeriod(); 

// initNecklace(1); // all ones 

// lastperiod = findSequencePeriod(); 


// debruijnsection[0] = 1; 

// debruijnsection[2*m + 1] = 0; 
// forGnti= 1;1<=m;i++) 

if { 

// — debruijnsection[i] = 0; 

i ft 

// forGnt i= m+1; i <= 2*m; i++) 
// { 

// — debruijnsection[i] = 1; 


} 
// } // end if 
// else 
// { 
for(int 1 = 0; 1 < p2; i++) 
{ 
debruijnsection[i + pl] = middlenecklace[i]; // middle third of deBruijn section 
} 
for(int i = 0; 1 < p3; i++) 
{ 
debruijnsection[i + pl + p2] = lastnecklace[i]; // last third of deBruijn section 
} 
// } // end else 
} // end createDeBruijnSection 
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// Method createWraparoundDeBrujjnSection 

// Return value none 

// Parameters none 

// Purpose concatenates nonperiodic portions of four adjacent necklaces: 
// 100, OOF O41 T1L21@ 
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void Necklaces::createWraparoundDeBruijnSection() 


{ 
wrapdebruijnsection[0] = 1; 
wrapdebruijnsection[2*m + 1] = 0; 


for(int i = 1; 1<=m; i++) 


{ 
wrapdebruijnsection[i] = 0; 
} 
for(Gint 1 = m+1; 1 <= 2*m; i++) 
{ 
wrapdebruijnsection[i] = 1; 
} 


} // end createWraparoundDeBrujnSection 
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// Method getInputDec 
// Return value none 
//Parameters long long input 
// Purpose obtains user input decimal value that yields binary string 
// to be searched for in DeBruijn sequence 
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void Necklaces::getInputDec(long long input) 
i 
inputstrdec = input; 
} //end getInputDec 
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// Method searchDeBrujjnString 
// Return value long long position 
//Parameters none 
// Purpose searches for user input binary string in DeBruijn sequence 
// and returns its location 
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long long Necklaces::searchDeBruijnString() 
{ 
long position = 0; // position from beginning of 111...1 necklace 
absposition = 0; 
int matchcount = 0; // counts how many matches there are between n-tuple and 
deBruijn string 
long stepcount = 0; // counts how many step you take thru deBruijn section until there 
is a match 
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cout << "\nsumperiod = '' << sumperiod << endl; 
for(int 1 = 0; 1 < sumperiod; i++) 
for(int j = 0; j <_m; j++) 


{ 


if(inputstring[j] == debruijnsection[i+j]) 


matchcount++; 

} 

else 

{ 
matchcount = 0; //reset 
break; 

} 


} // end inner for loop 


if(matchcount = = m) //section of deBruijn sequence fully and uniquely matches user 
input 
{ 
break; 
} 
stepcount++; // keep counting number of steps needed to find a match 


} 


position = stepcount + totaldistance; 


if(inputstrdec = = generatePower2(m)-1) // accounts for all ones necklace 


{ 


absposition = generatePower2(m); //last position in deBruin sequence read from right 
to left 


else if (inputstrdec = = 0) //accounts for all zeros necklace 

absposition = 1; //first position in deBruin sequence read from right to left 
ii //everything else in between two extremes 

absposition = generatePower2(m) - position; 


} 


return position; 
} // end searchDeBrujjnString 
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// Method findSequencePeriod 
// Return value int period 
//Parameters none 


// Purpose find period of a necklace 
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int Necklaces::findSequencePeriod() 


{ 


int period = 0; 
findDivisors(m); // need to find divisors first 


for(int 1 = 0; 1 < m; i++) 
{ 

if(divisors[i]!=0) 

{ 


for(int j = 0; j < (m/ divisors[i]) - 1; j++) 


period = 0; //reset 
for(int k = 0; k < divisors[i]; k++) 


{ 


if(sequence[k] == sequence[k + (j+1)*divisors[1i]]) 


period++; // period value keeps accumulating as long as there is a match 


; 


else 


{ 
period = m; 
break; // exit k for loop if there is no match 


} 
}// end k for loop 
if(period == m) 


break; // do not bother comparing remaining instances 


} 
} // end j for loop 


} // end outer if statement 
if(period == divisors[i]) 
break; // no need to check remaining divisors since we found period 
} 
} // end i for loop 


return period; 
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} // findSequencePeriod 
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// Method shiftSequenceToNecklace 

// Return value none 

//Parameters none 

// Purpose given an arbitrary n-tuple, it rotates the sequence until it generates 


necklace 
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void Necklaces::shiftSequenceToNecklace() 
{ 
int w = 0; // value of largest run of ones for cyclical rotation 
int wl = 0; // window equal to running ones without wraparound 
int w2 = 0; //first half of wraparound window 
int w3 = 0; // second half of wraparound window 


int firstsum = 0; // number of ones from i to m-1 when i+w-1 > m-1 
int secondsum = 0; // number of ones from 0 to i+w-1 

int onessum = 0; 

long long temp1 = 0; 

long long temp2 = 0; 


temp1 = binToDec(); // initial decimal value 
wl =countRunningOnes(0,m-1); //window size consisting of all ones 


for(int i = 0; 1< wl; i++) 
{ 

if(sequence[i]==1) 

{ 


w2++; 


} 


else 


{ 
break; 


} 
} 


for(int 1= m-1; 1 >= m-w]; i--) 
{ 
if(sequence[i]==1) 
{ 


w3+4+; 


} 
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{ 
break; 

} 
} 
if(wl > w2+w3) // need to find largest window size under cyclical rotation 
{ 

w=wl; 
} 
else 
{ 

Ww = w24+w3; 
} 


for(int 1 = 0; 1 < m; i++) 
{ 
if(i+w-1 <_m-1) // checks if end of window has reached end of sequence 
{ 
if(countRunningOnes(i,i+w-1)== w) 
{ 
decToBin(temp1); 
cycleshiftSequence(i); 
temp2 = binToDec(); // converts rotated sequence into a decimal value 
if(temp2 > temp1) 
{ 
temp1 = temp2; //stores only largest decimal equivalent associated with "w" 
ones 


} // end if(itw-1 <= m-1) statement 
else // may not need to use this 
{ 
continue; // ignores all windows which have less than "w" running ones 
} 
} 
else if (i+w-1 >= m-1) //window exceeds sequence length and wraps around to 
beginning 
{ 
//cout << "\nInside if('" << itw-1 <<" >=" << m-1 <<" ) statement..." << endl; 
//cout << "\n 1=" <<i<< endl; 
//cout <<'"\n temp! before cycleshift = " << temp1 << endl; 
//cout << "\n Binary value of temp1: "; 
decToBin(temp1); 
//printSequence(); 
//cout << "\n" << endl; 
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cycleshiftSequence(i); 
temp2 = binToDec(); // converts rotated sequence into a decimal value 
//cout << "\n Binary value of temp1 after shifting by "<<i<<":"; 
//printSequence(); 
//cout << "\n"; 
//cout << "\n temp2 =" << temp2 << endl; 
//cout << "\n temp! =" << temp! << endl; 
if(temp2 > temp1) 
{ 
//cout << "\n_ "; 
//cout << temp2 <<" >" << temp! << endl; 
temp1 = temp2; //stores only largest decimal equivalent associated with "w" ones 
} 


} 
}// end for loop 


decToBin(temp1); // converts decimal equivalent to binary 
} //end shiftSequenceToNecklace 
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// Method generateShiftingOnesSequence 

// Return value none 

// Parameters long long neckdec, decimal equivalent of necklace 


// long long numclump, number of missings after necklace 

// long distance, distance between adjacent signposts 

// long necknum, necklace number associated with each signpost 
// Purpose generates shifting ones sequence to use as signposts 


// Note: Entries are stored in a matrix format that is read from left 
// to right where the last entry of the previous row precedes 
// the first entry of the current row. This ensures I have enough 


// memory to collect enough signposts 
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void Necklaces::generateShiftingOnesSequence(long long necdec, long long clmpsz, 
long distance, long necknum) 
{ 

long long nval = necdec; 

long long nclump = clmpsz; 

long dist = distance; 

long nent = necknum; 


if(colindex <= MAX_COLUMNS) 
{ 


shiftonesarray[rowindex][colindex] = nval; 


107 


shiftonesclumpsizearray [rowindex][colindex] = clmpsz; 
shiftonesdistancearray[rowindex][colindex] = distance; 
shiftonesncountarray[rowindex][colindex] = ncnt; 
colindex++; // increment column index 
} 
else 
{ 
colindex = 0; // reset column index 
rowindex++; //increment row index 
if(rowindex <= MAX_LENGTH) 
{ 
shiftonesarray[rowindex][colindex] = nval; 
shiftonesclumpsizearray[rowindex][colindex] = clmpsz; 
shiftonesdistancearray[rowindex][colindex] = distance; 
shiftonesncountarray[rowindex ][colindex] = ncnt; 
colindex++; // increment column index 


} 


}// end generateShiftingOnesSequence 
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// Method getWeight 
// Return value void 
//Parameters int num 
// Purpose obtain user weight input for shifting ones sigposts 
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void Necklaces::getWeight(int num) 
{ 
weight = num; 
} //end getWeight 
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// Method getNumEvenlySpacedSignposts 
// Return value void 
//Parameters int num 
// Purpose obtain number of evenly spaced sigposts from user 
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void Necklaces::getNumEvenlySpacedSignposts(int num) 
{ 
numevensignposts = num; 
}// end getNumEvenlySpacedSignposts 
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// Method computeDensity 

// Return value int nclass, largest number of runnning zeros 
// Parameters int start, int finish 

// Purpose counts the total number of ones in the sequence 


// ** Note: To compute number of zeros = length of sequence - computeDensity 
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int Necklaces::computeDensity(int start, int finish) 


{ 
int alpha = start; 
int omega = finish; 
int density = 0; 


for(int 1 = alpha; 1 <= omega; 1++) 
{ 
if(sequence[i] == 1) 


{ 


density++; 
} 
} 


return density; 
} //end computeDensity 
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// Method testforNecklace 

// Return value bool foundnecklace 

// Parameters long long bininput 

// Purpose tests n-tuple to see if it is a necklace 
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bool Necklaces::testforNecklace(long long bininput); 

{ 

bool foundnecklace = false; 

long long ndecimal = bininput; // binary decimal associated with n-tuple 

nstring = decToBin(ndecimal); // subroutine that converts decimal value into a binary 

// sequence 

shiftSequenceToNecklace(); // subroutine that shifts a binary sequence into its necklace 

// representative 


shiftdecimal = binToDec(); // subroutine that converts a binary string into its decimal 
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// equivalent 


if(shiftdecimal = = ndecimal) // test to see if current n-tuple is a necklace 


{ 

foundnecklace = true; 
} 
else 
{ 

foundnecklace = false; 
} 


return foundnecklace; 


// end testforNecklace 
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// Method countNecklacesPerClass 

// Return value none 

//Parameters int ones, number of leading ones 

// Purpose separates the number of leading ones in each necklace 


// into distinct classes 
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void Necklaces::countNecklacesPerClass(int ones) 


{ 


int p = ones; //number of leading ones 
if(p != m-q)//checks if current ones count is different from last ones count 


onescount[m-p]= 1; // ensures unique number of ones is counted only once 
q++; 
} 


else 


{ 


onescount[m-p]++; //increase number of necklaces with same number of leading ones 


} 


if(p == 0) //accounts for all zeros string 


{ 


onescount[m-p]= 1; 


} 


}// end countNecklacesPerClass 
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// Method generateChangeSequence 

// Return value none 

//Parameters long long oldec, decimal value of previous necklace 


// long long newdec, decimal value of current necklace 

// int indx, array index for sequence 

// Purpose counts the number of bits that differ between each necklace 
// and inserts these values into a sequence 
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void Necklaces::generateChangeSequence(long long olddec, long long newdec, int indx) 
{ 

int oldseqjMAX_LENGTH]; // binary sequence for previous necklace 

int newseqi[MAX_LENGTH]; // binary sequence for current necklace 

int bitchange = 0; // counts the number of bit changes between both sequences 

int idx = indx; // column index for bit change sequence 


decToBin(olddec); 
for(int i=0; 1 < m; i++) 
{ 
oldseq[i] = sequence[i]; 
} 
decToBin(newdec); 
for(int i=0; 1 < m; i++) 
{ 
newseq[i] = sequence[i]; 


} 
for(int i=0; 1 < m; i++) 


if(oldseq[i]!= newseq[1]) 
{ 
bitchange++; 
} 
} 
bitchangeseq[idx]= bitchange; 
}// end generateChangeSequence 


1/5 7 2 2 2 2g 2g 2g ae ae ae he he he he fe he he fe he he 26 26 26 28 28 28 2s 2s 2s 2s 2 2 2 2 2 2 28 2 28 
// Method binToDec 

// Return value int dec 

// Parameters none 


// Purpose converts binary number of a sequence into its decimal equivalent 
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long long Necklaces::binToDec() 


long long dec = 0; 
for(int 1 = 0; 1< m; i++) 
{ 
dec = dec + generatePower2(m- | -i)*sequence[i]; 
} 
return dec; 
}// end binToDec 
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// Method decToBin 
// Return value none 
//Parameters _long long decimal input 


// Purpose converts decimal input into its binary equivalent sequence 
[7 7 2 2 2g 2g 2g ae ae ae ae ee he he fe he fe he he ie ie ie 26 28 28 26 2s 2s 2s 2s 2s 2 2 2 2 2 2 2 2k 


void Necklaces::decToBin(long long input) 
{ 
int r = 0; //number of leading ones 
int j = m-1; //start at LSB 
initNecklace(r); 
long long quotient = input; //find a way to change to long long 
while(quotient !=0 && j >= 0) 


{ 
if (quotient == () 
{ 
sequence[j] = 1; // accounts for the MSB where 1/2 = 0 + (rem=1) 
} 
else 
{ 


sequence[j] = quotient%2; //remainder should be 0 or 1 


} 
J-> 
quotient = quotient /2; 


} 
}// end decToBin 
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// Method createUnobservedSequences 

// Return value none 

// Parameters int decl, int dec2 

// Purpose generates preabsorbed necklaces that are skipped by theta 
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// and which exist between the last necklace in a given class 


// and the first necklace in the next adjacent class 
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long Necklaces::createUnobservedSequences(int decl, int dec2) 

{ 
long index = 0; 
long index1 = dec1; // decimal equivalent of last necklace in a given class 
long index2 = dec2; // decimal equivalent of first necklace in next class 
long seqcount = 0; //counts sequences 


//index = index1-2; 

index = index]; 

while (index >= index2) 

{ 
decToBin(index); //generates binary equivalent of decimal 
printSequence(); 
cout << "decimal value = '' << index << endl; 
index = index - 2; // substract 2 to skip over binary strings ending in | 
seqcount++; 
cout << "\n"; 

} 

return seqcount; 

} // end createUnobservedSequences 
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// Method findDivisors 
// Return value none 
//Parameters int num 


// Purpose find all divisors for an integer n 
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void Necklaces::findDivisors(int num) 


{ 
for(int i=0; i< MAX_LENGTH+1; i++) 
{ 
divisors[i] = 0; 
} 


int number = num; 
for(int i=0; i < number; i++) // does not include m itself 


{ 


if(Mumber%(i+1)== 0) // checks for divisors; i+1 ensures no division by zero 


{ 
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divisors[i] = i+1; // assigns divisors to the (1-1)th place in array 
} 


} 
}// end findDivisors 


J [7B AR AR A A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee he he fe he fe he he ie Ie ie ie 2s Ie iS 


// Method printDivisors 
// Return value none 
// Parameters int num 


// Purpose prints all divisors for an integer n 
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void Necklaces::printDivisors(int num) 

{ 
int number = num; 
cout << '\nDivisors of '' << number << "...\n'' << endl; 
for(int i=0; 1 < number; i++) 


{ 
if(divisors[i]!=0) 
{ 
cout << setw(3) << divisors[i] << " "'; // divisors printed 
} 
} 


} // end printDivisors 


J [FR FR A AR A A 2 2 2 2g 2 2 2 2 2 2 2 2 2 2 2 2 3 2 ae ee he fe he fe he he ie Ie ie ie is Ie is 


// Method findGCD 
// Return value int gcd 
//Parameters int numl, num2 


// Purpose find greatest common divisor (gcd) for two integers 
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int Necklaces::findGCD(int num1, int num2) 
{ 

int first = num 1; 

int second = num?; 

int factorl1[MAX_ LENGTH + 1]; 

int factor2,;/MAX_ LENGTH + 1]; 


for(int i=0; i< MAX LENGTH + 1; i++) 
{ 

factor] [i] = 0; 

factor2[i] = 0 
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} 


int factorlcnt = 0; 

findDivisors(first); 

for(int i=0; i< MAX LENGTH + 1; i++) 
{ 


if(divisors[i]!=0) //only collect non-zero divisors 


factor1[i]= divisors[i]; 
factorlcnt++; 
} 
} 


int factor2cnt = 0; 
findDivisors(second); 
for(int i=0; i< MAX LENGTH + 1; i++) 
{ 
if(divisors[i] !=0) //only collect non-zero divisors 
{ 
factor2[i]= divisors[i]; 
factor2cnt++; 
} 
} 


int factorcnt = 0; 
if(factorlcnt > factor2cnt) 


{ 

factorcnt = factorl cnt; 
} 
else 
{ 

factorcnt = factor2cnt; 
} 
int gcd = 0; 


for(int 1=0; 1 < factorcnt; i++) 
for(int j=0; j < factorcnt; j++) 
if(factor2[j] == factor 1 [i]) 
if(factor2[j] > gcd) 
ged = factor2[j]; 
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} 
} 
} 
} 


return gcd; 
} // end findGCD 
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// Method eulerTotient 
// Return value long long Zn 
//Parameters int num 


// Purpose Calculate number of necklaces for a given n 
J [7 FR A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee fe he fe fe fe he he ie ie ie ie Ie ie 


long long Necklaces::eulerTotient(int num) 
//double Necklaces::eulerTotient(int num) 
{ 

int size = num; 

long long Zn = 0; 

long long totientiMAX LENGTH+1]; // contains number of integers relatively prime 
to the divisor including "1" 

//double Zn = 0; 

//double totientt MAX _LENGTH+1]; // contains number of integers relatively prime to 
the divisor including "1" 


int divisorarray[MAX_LENGTH+1]; 


for(int i=0; 1<= MAX_LENGTH; i++) 
{ 

totient[i] = 0; 

divisorarray[1i] = 0; 


} 
int divisorcnt = 0; 
findDivisors(size); 
for(int i=0; i <= MAX_LENGTH; i++) 
{ 
if(divisors[i]!=0) //only collect non-zero divisors 
{ 
divisorarray[i] = divisors[i]; 
divisorcnt++; 
} 
} 
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for(int i=0; 1 <= size; i++) 


if(divisorarray [i] != 0) 


{ 
for(int j = 1; j < divisorarray[i]; j++) 
{ 
if(findGCD(divisorarray[i],j)== 1) 
{ 
totient[i]++; 
} 
} 
} 
} 


totient[O] = 1; // to include "1" 
long long sum = 0; 
//double sum = 0; 
for(int i=0; 1 <= size; i++) 
{ 
if(totient[1]!=0) 
{ 
sum = sum + totient[i]*generatePower2(m/divisorarray[1]); 
} 
} 


Zn = sum / size; 


cout << '\n"; 
return Zn; 
}// end eulerTotient 
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// Method generatePower2 
// Return value long long power2val 
//Parameters long long exponent 


// Purpose generates 2 raised to any integer 
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long long Necklaces::generatePower2(long long exponent) 


{ 


long long power2val = 1; 
long long expint = exponent; 
for(int 1 = 0; 1 < expint; 1++) 
{ 


power2val *= 2; //calculates power of 2 
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} 


return power2val; 


}// end generatePower2 
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// Method displayNecklaceInfo 


// Return int ncountarray 
//Parameters none 
// Purpose displays all information related to a necklace 
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void Necklaces::displayNecklaceInfo(int a, int b, long long c, long long d) 
{ 

int column = a; 

int clssnum = b; 

long long val = c; 

long long totalclumps = d; 


cout << '\nclass number = " << clssnum << end]; 

cout << '\nnumber of clumps so far = " << totalclumps << endl]; 

cout << "\nclumpcentarray[" << clssnum << '']= '"' << totalclumps << endl; 

cout << ''\nncountarray[''<< clssnum << "]["' << column << "] =" << ncount << 
endl; 

//cout << "\ndecimalarray["<< clssnum << "][" << column << "] =" << 
decimalarray[clssnum][column] << endl; 

cout << '\n"; 
}// end displayNecklaceInfo 
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// Method displayPower2Array 

// Return value none 

//Parameters none 

// Purpose displays power2[i] array 
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void Necklaces::displayPower2Array() 


{ 
int al = 3; 
int a2 = 2; 
for(int i=al; i<=m; i++) 
{ 
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for(int j=a2; j<1; j++) 


power2[i] *= 2; //calculates power of 2 


} 
} 


cout << '"\n\nPower of 2 approximation for each class number:\n"' << end]; 
for(int i=0; i <= m; i++) 
{ 

cout << setw(OUTPUT_CELL_SIZE)<< power2[i] <<"""'; 


} 
}// end displayPower2Array 
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// Method displayDifference 

// Return value none 

//Parameters none 

// Purpose calculates and displays difference between power2 approximation and 


number of necklaces 
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void Necklaces::displayDifference() 
{ 


for(int i = 0; 1 <= m; i++) 


diff[i] = power2[i] - onescount[i]; 


} 


cout << '\n\nNumber of differing necklaces per class number:\n" << end]; 
for(int i=0; i <= m; i++) 
{ 

cout << setw(OUTPUT_CELL_SIZE)<< diff[i] <<""'"'; 


} 
}// end displayDifference 


J [RFR A i 2 2 2 2 2 2 He 2 2 2 2 2 2 2 2h 2 2h 2 2 2 2 he fe he fe he fe he ie ie ie ie ie Ie is 
// Method display NumbNeck 

// Return value none 

// Parameters none 


// Purpose displays number of necklaces per class 
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void Necklaces::display NumbNeck() 


119 


{ 
cout << '\nNumber of necklaces per class number:\n" << endl; 
for(int i=0; 1 <= m; i++) 
{ 
cout << setw(OUTPUT_CELL_SIZE)<< onescount[i] <<" "'; 


} 
}// end displayNumbNeck 
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// Method displayClassNumbers 
// Return value none 
// Parameters none 


// Purpose displays class numbers of necklaces in decreasing order 
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void Necklaces::displayClassNumbers() 
{ 
cout << '\nDecreasing class numbers:\n" << end]; 
for(int 1 = m; 1 >= 0; i--) 
{ 
cout << setw(OUTPUT_CELL_SIZE) <<i<<"'"; // displays decreasing leading ones 
} 
cout << endl; 
}// end displayClassNumbers 
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// Method displayNcount 
// Return value none 
// Parameters none 


// Purpose displays total number of necklaces 
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void Necklaces::displayNcount() 
{ 
cout << '\n\n*** Total number of theta generated binary " << m << '"'-tuples = '"' 
<< ncount + noncount << '"' ***" <<endl; 
cout << ''\nNumber of necklaces = '' << ncount << endl; 
cout << '\nNumber of theta generated non-necklaces = '' << noncount << endl; 
cout << '"\nNumber of missing even '' << m << ''-tuples = "' << totmissing << endl; 
cout << '"\nNumber of missing even '' << m << ''-tuples + number of necklaces = " 
<< totmissing + ncount << endl; 
cout << '\nMaximum distance of necklaces = "' << distance << endl; 
cout << '"\n"; 
}// end displayNcount 
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// Method displayClumpsPerClass 
// Return value none 
// Parameters none 


// Purpose displays total number of missing even n-tuples per class 
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void Necklaces::displayClumpsPerClass() 
{ 
long clumptotal = 0; 
cout << '\n\nNumber of missing clumps per class number:\n" << end]; 
for(int 1 = m; 1 >= 0; i--) 
{ 
cout << setw(OUTPUT_CELL_SIZE)<< clumpcentarray[i] <<" "; 
clumptotal += clumpcntarray[i]; 
} 
cout << '\n\nTotal number of clumps = " << clumptotal << "' for threshold = '' << 
threshold << endl; 


}// end displayClumpsPerClass 
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// Method displayClumpDistribution 
// Return value none 
// Parameters none 


// Purpose displays various clump sizes per class 
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void Necklaces::displayClumpDistribution() 
t 

cout << "\n\n\nDistribution of missing necklaces for n = " << m <<" and threshold 
= '' << threshold << ":'' << endl; 

cout << '"\nEach number below represents the size of the clump." << end]; 

cout << '\n"; 


long long summissing = 0; 

long long nummiss = 0; 

int numclump = 0; 

int listwidth = 4; // displays only so many triples across 

// triple = (Necklace Number, Decimal Equivalent of Necklace, Number of Missings 
after Necklace) 


// percent = percentage of the maximum clump size 
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cout << "Clump sizes which are greater than " << num << " percent of the 
maximum clump size are shown. " << end]; 

cout << "\nMaximum clumpsize for all the classes = '' << maxclumpsize << endl]; 

cout << '\nNote: Lower necklace decimal value = Upper necklace decimal value - 
2*(number of missings)\n"' << endl; 


cout << '\n"; 
for (int i = 0; 1 < m; i++) 
{ 


cout << "\nClass '"<<i<<'"': "; 
summissing = 0; //resetting count 
nummiss = 0; 

for (int j = 0; j < MAX_COLUMNS; j++) 


if(clumpsizearray[1i][j]!= 0) 
{ 
if(clumpsizearray[1i][j] >= maxclumpsizearray[i]*percent) 
{ 
cout << setw(OUTPUT_CELL_SIZE) << clumpsizearray[i][j] << "'"'; 
summissing += clumpsizearray[1][j]; 
nummiss ++; 
} 
} 
} 
if(nummiss!=0) 
{ 
cout << '\n\nTotal number of missing even " << m << "'-tuples exceeding 
threshold of " << threshold << '' = " << summissing << endl; 
cout << '"\nTotal number of missing even '' << m << ''-tuples below threshold 
of '' << threshold << '' = '' << missingbelowthreshold[i] << endl; 
cout << '\nTotal number of clumps = " << nummiss << end]; 
cout << "\nMaximum clump size for this class = '' << maxclumpsizearray[i] << 
endl; 
cout << '\n(Necklace Number, Decimal Equivalent of Necklace, Number of 
Missings after Necklace)''<< endl; 
cout << '\n"; 
numclump = 0; //resetting count 
for (int j = 0; }] < MAX _COLUMNS; j++) 


/Af(ncountarray[i][j] != 0 && missingarray[i][j]!= 0 && decimalarray[i][j]!= 9) 
{ 
//cout << setw(OUTPUT_CELL_SIZE) << "(" << ncountarray[i][j] <<" ," << 
decimalarray[i][j]<< " , " << missingarray[i][j] << ")" ; 
numclump++; 
if(j+1)%listwidth == 0) 
{ 
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cout << "\n\n"; 


} 
} 

} // end for 
} // end if 
else 
{ 

cout << ''No missings in this class.'' << endl; 
} 
cout << "\n\n" << endl; 

}//end outside for 
}// end displayClumpDistribution 
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// Method displayDecimalArray 
// Return value none 
// Parameters none 


// Purpose displays decimal equivalent of necklaces associated with clumps per class 
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void Necklaces::displayDecimalArray() 
{ 

int listwidth = 4; // displays only so many triples across 

cout << "\n\n\nDistribution of decimal equivalents for missing clumps for n =" << 
m << '" and threshold = "' << threshold << '':" << endl; 

cout << '"\n"; 

for (int i = 0; 1 < m; i++) 


for (int j = 0; j < MAX_COLUMNS; j++) 


{ 
/Af(decimalarray[i][j]!= 9) 
{ 
//cout << setw(OUTPUT_CELL_SIZE) << decimalarray[i][j]<<"" ; 
if((j+1)%listwidth == (0) 
{ 
cout << '"\n\n"; 
} 
} 
} 
} // end for 


cout << "\n\n\nList of necklaces associated with clumps for n = '' << m <<" and 
threshold = '' << threshold << ":' << endl; 
cout << '\n"; 
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for (int i= m-1; 1 >= 0; i--) 


for (int j = 0; j < MAX_COLUMNS; j++) 


{ 
/Af(decimalarray[i][j]!= 9) 
{ 
//decToBin(decimalarray [i][j]); 
printSequence(); 
cout << '\n\n"; 
} 
} 
} 


}// end displayDecimalArray 
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// Method displayShiftingOnesSignposts 

// Return value none 

// Parameters none 

// Purpose displays decimal and binary equivalents of shifting ones signposts 


// along with clump sizes 
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void Necklaces::displayShiftingOnesSignposts() 
{ 
int numshiftones = 0; // total number of shifting ones signposts 
int loopcnt = 0; // keeps track of long long for loop count 
long long shiftonesmissingsum = 0; //sums missings 
char answerl1; 
char answer2; 
cout << ''\n\n\nSignposts with shifting ones for n = '' << m << endl]; 
cout << "\n"; 


cout << 'Do you want to see the signposts? (Enter 'y' or 'n'): "; 
cin >> answerl; 
cout << '\n"; 


cout << 'Do you want to see the various clump sizes? (Enter 'y' or 'n'): "; 
cin >> answer?2; 
cout << '"\n"; 


for(int 1 = m-1; 1 >= 0; 1--) 

{ 
for (int j = 0; |] < MAX_COLUMNS; j++) 
{ 
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if (shiftonesarray[1i][j]!= mark) 


{ 
numshiftones++; 
if(answerl == 'Y' || answerl == 'y') // displays binary strings, decimal values and 
clump sizes 
{ 
decToBin(shiftonesarray [1][j]); 
printSequence(); 


cout <<" decimal value = "<< shiftonesarray[i][j] <<", "5 
cout << setw(OUTPUT_CELL_SIZE) << "clump size = ''<< 
shiftonesclumpsizearray[i][j] << endl; 
shiftonesmissingsum += shiftonesclumpsizearray[1][j]; 
cout << '\n"; 
} 
else // displays only various clump sizes 
{ 
loopent++; 
if (shiftonesclumpsizearray [i][j]!= 0) 
{ 


if(answer2 == 'Y' || answer2 == 'y') 


{ 
if doopcnt <= 1) // ensures output statement is not repeated more than once 
{ 
cout << "Clump sizes for shifting ones signposts: '' << endl; 
cout <<'"'\n"'; 
} 
cout << setw(OUTPUT_CELL_SIZE) << shiftonesclumpsizearray[i][j] <<" 
} 
shiftonesmissingsum += shiftonesclumpsizearray[1]|[j]; 
} 
} 
} 
else 
{ 
continue; 
} 
} 


} 

binlength = generatePower2(m- 1); 
int k1 = 4; // equipment cost 

int k4 = 1; // temporal cost 


petriumilestones = floor(binlength / m); // assume max period = m 
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milestonedistance = floor(binlength /shiftonesclump); // distance between milestones 
using number of shifting ones signposts 
topt = ceil(sqrt((k1/k4)*binlength));// gives bad answer 


cout << '\n\nNumber of shifting ones signposts associated with missings (using 
array) = '' << numshiftones << endl; 

cout << '\nNumber of shifting ones signposts associated with missings (using pure 
count) = '' << shiftonesclump << endl; 

cout << '\nNumber of shifting ones signposts NOT associated with missings = '' << 
shiftonesnoclump << endl; 

cout << "\nTotal number of missings accounted for by shifting ones = '' << 
shiftonesmissingsum << end]; 

cout << "\n\nPetriu length is 24"' << m <<" - 1 =" << generatePower2(m) - | << 
endl; 

cout << '\nTotal number of Petriu's signposts (assuming the max separation 
distance is '' << m <<"' bits) = '' << petriumilestones << endl]; 

cout << ''\nPetriu distance between milestones (assuming the number of milestones 
is '' << shiftonesclump << '') = "' << milestonedistance << endl; 

cout << "\nOptimal Petriu distance between milestones = '' << topt << endl; 
}// end displayShiftingOnesSignposts 
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// Method displayShiftingOnesDistances 
// Return value none 
// Parameters none 


// Purpose displays distances between shifting ones signposts 
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void Necklaces::displayShiftingOnesDistances() 
{ 
long reldistance = 0; // relative distance 
long currentdistance = 0; 
long sumdistance = 0; // sum of rleative distance elements 
long numdistance = 0; // number of distance elements 
long avedistance = 0); 
int colspread = 15; // width of data elements 
int spreadcount = 0; // number of data elements used to spread data matrix evenly 
long shiftonespetriunumbigger = 0; // number of distances greater than petriu distance 
long shiftonespetriunumlesser = 0; // number of distances less than or equal to petriu 
distance 


cout << " "'; // adjustment for first row to make entries look even 
for(int i = 0; 1< MAX_LENGTH,;; i++) 
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{ 
for(int j = 0; } < MAX_COLUMNS; j++) 


if(shiftonesdistancearray[i][j] != 0) 

{ 
reldistance = shiftonesdistancearray[i][j] - currentdistance; 
cout << setw(OUTPUT_CELL_SIZE) << reldistance; 
currentdistance = shiftonesdistancearray[i][j]; 
sumdistance += reldistance; 
numdistance++; 
spreadcount++; 
//if(reldistance > milestonedistance) 
if(reldistance > topt) 


{ 
shiftonespetriunumbigger++; 
} 
else 
{ 
shiftonespetriunumlesser++; 
} 
if(spreadcount == colspread) 
{ 
cout << '"\n"; 
spreadcount = 0; //reset 
} 
else 
{ 
cout << setw(OUTPUT_CELL_SIZE) << ''"'; 
} 


} 
} 
} 


avedistance = sumdistance / numdistance; 
cout << '\n\navedistance = "' << avedistance; 
//cout << "\n\nNumber of signposts whose distance is greater than " << 


milestonedistance << " =" << shiftonespetriunumbigger << endl; 
//cout << "\n\nNumber of signposts whose distance is less than or equal to " << 
milestonedistance << " =" << shiftonespetriunumlesser << end]; 


cout << "\n\nNumber of signposts whose distance is greater than "' << topt <<" =" 
<< shiftonespetriunumbigger << endl; 

cout << '\n\nNumber of signposts whose distance is less than or equal to " << 
topt<< " = '' << shiftonespetriunumlesser << endl; 


}// end displayShiftingOnesDistances 
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// Method displayBitChangeSequence 
// Return value none 
// Parameters none 


// Purpose displays bit changes between necklaces 
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void Necklaces::display BitChangeSequence() 
{ 


cout << ''\n\nThe following sequence shows the number of bit changes between 
each pair of necklaces. '"' << endl; 
cout << '"\n"; 


for(int i=0; 1< MAX_COLUMNS; i++) 
{ 
if (bitchangeseq[i]!= 0) 


cout << setw(OUTPUT_CELL_SIZE) << bitchangeseq|[i] << '"""; 
} 


} 
}// end displayBitChangeSequence 
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// Method displayFirstNecklace 
// Return value none 
// Parameters none 


// Purpose displays prior necklace in deBruijn section 
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void Necklaces::displayFirstNecklace() 
for(int 1 = 0; 1 < m; i++) 


if(firstnecklace[i]!= 9) 
{ 


cout << firstnecklace[i] << '"' "; 


} 
} 
cout <<"\n"'; 
cout << '\nfirstperiod = '' << firstperiod << end]; 


}// end displayFirstNecklace 
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// Method displayMiddleNecklace 
// Return value none 
// Parameters none 


// Purpose displays input necklace in deBruijn section 
J [FRR A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee fe he fe he he he ie ie ie ie ie 2S it 


void Necklaces::displayMiddleNecklace() 
{ 
for(int 1 = 0; 1 < m; i++) 
{ 
if(middlenecklace[i]!= 9 ) 
{ 
cout << middlenecklace[i] << 
} 
} 


cout <<'\n"; 
cout << "\nmiddleperiod = '' << middleperiod<< endl; 


rom, 


}// end displayMiddleNecklace 
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// Method displayLastNecklace 
// Return value none 
// Parameters none 


// Purpose displays post necklace in deBruijn section 
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void Necklaces::displayLastNecklace() 
{ 
for(int 1 = 0; 1< m; i++) 
{ 
if(lastnecklace[i]!= 9 ) 
{ 
cout << lastnecklace[i] << 
} 
} 


cout <<"\n"'; 
cout << '\nlastperiod = '' << lastperiod << endl; 
cout << "\nsumperiod = '' << sumperiod << endl; 


wom, 


}// end displayLastNecklace 
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// Method displaydeBruijnSection 

// Return value none 

// Parameters none 


// Purpose displays deBruijn sectional associated with user input necklace 
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void Necklaces::displaydeBruijnSection() 


for(int 1 = 0; 1 < firstperiod; 1++) 
if(debruijnsection[i]!= 9) 
cout << debruijnsection [i] << '' ''; // prints the first third 
} i end for loop 
cout << "4 '"'; 


for(int 1 = 0 ; 1< middleperiod ; i++) 
{ 


if(debruijnsection[i + firstperiod]!= 9) 


{ 


cout << debruijnsection [i + firstperiod] <<" ''; // prints the second third 


i 
} // end for loop 
cout << "4 '"'; 


for(int 1 = 0; 1 < lastperiod; i++) 
if(debruijnsection[i + firstperiod + middleperiod]!= 9) 


{ 


cout << debruijnsection [i + firstperiod + middleperiod] << " "'; // prints the last 
third 


} 
} // end for loop 
//\ // end else 


cout << "\n''<< endl; 


} // end displaydeBruijnSection 
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// Method display UnobservedSequences 
// Return value none 
// Parameters none 
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// Purpose displays all missings (even binary n-tuples) between any 2 nunmbers 
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void Necklaces::display UnobservedSequences() 
{ 
char answer; 
int input] =m ; //need to change to long long 
int input2 = 0; //need to change to long long 
long totseq = 0; 


cout << '\n\nWould you like to see all the binary sequences between 2 numbers? (y 
orn)"; 
cin >> answer; 
if (answer == 'Y' |] answer == 'y') 
{ 
cout << ''\nThis section will help you see all the sequences between 2 numbers." 
<< endl; 
cout << ''\nPlease enter 0 to quit or an even positive number less than " << 
generatePower2(input1)-1 << '"': "; 
cin >> input1; 


while(input1%2 !=0) 
{ 
cout << ''\nPlease enter another positive number: "'; 
cin >> input; 
} 
if(inputl != 0) 
cout << ''\nPlease enter another number less than or equal to " << input! <<": "; 
cin >> input2; 
cout << ''\nGenerating unobserved sequences..."" << end]; 
cout << "\n"'; 
totseq = createUnobservedSequences(input1 ,input2); 
cout << ''\n\nTotal number of sequences = "' << totseq; 


} 


} 
}// end displayUnobservedSequences 
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// Method outputData 

// Return value none 

// Parameters none 

// Purpose provide numerical data useful for analysis 
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void Necklaces::outputData() 


{ 
displayNcount(); 


displayClassNumbers(); 
displayNumbNeck(); 
displayPower2Array(); 
displayDifference(); 
displayClumpsPerClass(); 
displayDecimalArray(); 
displayClumpDistribution(); 
displayShiftingOnesSignposts(); 
display BitChangeSequence(); 
display UnobservedSequences(); 
}// end outputData 


// end of file Necklaces.cpp 


2. Thesismain CPP File 
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// Thesismain.cpp 
// 
// LT John Ortiz 
// 
// 
// Project: Thesis 
// Operating Environment: Windows XP Home 
// Compiler: Visual Studio .NET 
// Date: 08 December 2005 
// Description: This program generates necklaces and performs analysis on them 
// 
// Inputs: This program takes as inputs the length of the necklace 
// Outputs: This program displays a table of necklaces 
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// Processes: none 

// Assumptions: none 
// 

// Warnings: none 

// 
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#include "Necklaces.h" 

#include "[OThesis.h" 

#include "SieveSizeException.h" 
#include "math.h" 

#include <iostream> 

#include <iomanip> 

#include <cctype> 

#include <fstream> 


using std::ostream; 
using std::ofstream; //output file stream 
using namespace std; 


J [FR PH i 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2h 2 2 2 2 he fe fe fe fe he he ie ie ie ie ie Ie it 


// Method checkInput(input2); 
// Return value int input2 
//Parameters int input2 
// Purpose checks if input2 is valid 
J [78 78 A AB A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ee he fe he he fe he fe he ie ie ie 2k 28 ik 
int checkInput(int input1, int input2) 
{ 
int oneslength = input2; 
int limit = input]; 
while (oneslength < 0 Il oneslength > limit) 
{ 
cout << "Please enter a number between 0 and " << limit << " for the initial 
sequence." << endl; 
cin >> oneslength; 
} 
return oneslength; 
} // end checkInput 


int main() 

I 
long long input1 = 0; 
long long input2 = 0; 
long long input3 = 0; 
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char input; 

char answer; 

long long inputstep = 0; 
int maxinput = 63; 


int num = 0; 
int numl = 0; 
int num2 = 0; 


int intshift = 0; 

int onesrun = 0; 

int decimal = 0; 

int necknum = 0; 

int cellsize = 3; 

int bound = 0; // value of threshold 

int length = 0; //length of binary string 
int permaxclump = 0; // percentage of the maximum clump size 
long long expinput = 0; 

long long expoutput = 0; 

bool flag = false; 


cout << "This program will help you generate a list of necklaces" << endl; 
cout << "to enable you to perform analysis on them." << end]; 
flag = false; 
while (flag == false) 
{ 
cout << "\n\nPlease enter an integer for your necklace length between | and 50: "; 
try 
{ 
length = IOThesis::inputInt(); // verifies user input dimension is valid 
flag = true; 
} 
catch(SieveSizeException &sieveSizeException) 
{ 
cout << "Error occured: " << sieveSizeException.what() << endl; 
flag = false; 
} 
} 


cout << "\n\nControlling the threshold (smallest size of missings displayed) and " << 
endl; 

cout << "\nthe percentage of the maximum clump size can help manage the limitations 
" << endl; 

cout << "\nof memory when displaying the missings. If you don't want to adjust the " 
<< endl; 

cout << "\nthreshold or percentage, enter '0'. Else, hit any key: "; 

cin >> answer; 
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if (answer != '0') 

{ 
cout << "\n\nPlease enter the threshold: "; 
cin >> bound; 


cout << "\nPlease enter the percentage: "; 
cin >> permaxclump; 


} 


Necklaces necklace = Necklaces(length, bound, permaxclump); //Necklace class 
instantiated 


// Relatively prime test 
cout << "\n" << endl; 
cout << "Would you like to know whether two integers are relatively prime? "; 
cin >> answer; 
while (answer == 'y' Il answer == 'Y') 
{ 
cout << "\nPlease enter your first integer: "; 
cin >> input1; 
cout << "\nPlease enter your second integer: "; 
cin >> input2; 
cout <<"\nThe greatest common divisor between " << inputl <<" and " << input2 << 
"is " << necklace.findGCD(input1 ,input2); 
if(mecklace.findGCD(input1 ,input2)==1) 
{ 
cout << "\n\nThe integers are relatively prime." << endl; 
} 
else 
{ 
cout << "\n\nThe integers are not relatively prime." << endl; 
} 
cout << "\n\nWould you like to know whether another two integers are relatively 
prime? "; 
cin >> answer; 
} // end relatively prime test 


//Calculating divsiors of binary string length 
cout << "\n" << endl; 


cout << "Would you like to know number of necklaces for and the divisors of n=" << 
length <<" ?"; 

cin >> answer; 

if (answer == 'y' || answer == 'Y') 

{ 


necklace.findDivisors(length); 
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necklace.printDivisors(length); 
cout << "\nCalculated value of number of necklaces = "<< 
necklace.eulerTotient(length) << endl; 


} 


cout << "\nWould you like run the Necklace algorithm? "; 
cin >> answer; 

if (answer == 'y' II answer == "Y') 

{ 

cout << "\nWould you like to generate signposts? "; 

cin >> answer; 


if (answer == 'y' || answer == 'Y') 
{ 
cout << "\n\nThis next section enables you to generate signposts using the shifting ones 
strategy."<< endl; 
cout << "\nPlease enter a 'l' if you want to use the shifting ones signposts: " ; 
cin >> answer; 
if (answer == '1') 
{ 
cout << "\n\nPlease enter the weight you would want to use: "; 
cin >> num; 
necklace.getWeight(num); 
} // end if 
cout << "\n" << endl; 


cout << "If you want to generate all the necklaces, enter '1' or else hit any other key: "; 
cin >> answer; 
if (answer !='l') 
{ 
cout << "\nPlease enter a number between 0 and " << 
necklace.generatePower2(length)-1 <<" for the intial decimal: "; 
cin >> input1; 
cout << "\nPlease enter a number between 0 and " << 
necklace.generatePower2(length)-1 <<" for the final decimal: "; 
cin >> input2; 
cout << "\nPlease enter a number between 0 and " << 
necklace.generatePower2(length)-1 <<" for the input decimal: "; 
cin >> input3; 


else //default settings to generate all the necklaces 


{ 
inputl = necklace.generatePower2(length)-1; 
input2 = 0; 
input3 = 0; 
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} 
cout << "\n\nPlease enter 'y' or 'Y' if you would like to see the necklaces. Otherwise, 
enter any key: "; 


cin >> input; 


if (input == 'y' Il input == "Y") 


{ 
flag = true; 
cout << "\nGenerating Necklaces..." << endl; 
} 
else 
flag = false; 
} 


necklace.runNecklaceAlgorithm(input1, input2, input3, flag); 
necklace.outputData(); 


} 


cout << "\n\nDo you want to test cycle shifting a sequence? "; 
cin >> answer; 


while (answer == 'Y' || answer == 'y’) 
{ 
cout << "\nPlease enter a decimal to generate a binary sequence: "; 
cin >> input1; 
necklace.decToBin(input1); 
cout << "\n"; 
necklace.printSequence(); 
cout << "\n"; 
cout << "\nPlease enter an integer to test the cyclic shift: "; 
cin >> intshift; 
necklace.cycleshiftSequence(intshift); 
cout << "\n"; 
necklace.printSequence(); 
cout << "\n"; 
cout << "\n\nDo you want to test cycle shifting another sequence? "; 
cin >> answer; 


} 


cout << "\n\nDo you want to compute the density of a sequence? "; 
cin >> answer; 
while (answer == 'Y' |I answer == 'y’) 


{ 


cout << "\nPlease enter a decimal to generate a binary sequence: "; 
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cin >> input1; 
necklace.decToBin(input1); 
cout << "\n"; 
necklace.printSequence(); 


cout << "\n"; 
cout << "Density =" << necklace.computeDensity(0, length-1) << endl; 


cout << "\n\nDo you want to compute the density of another sequence? "; 
cin >> answer; 


} 


cout << "\n\nDo you want to find the period of a sequence? "; 
cin >> answer; 
while (answer == 'Y' Il answer == 'y’) 


{ 


cout << "\nPlease enter a decimal to generate a binary sequence: "; 
cin >> input1; 

necklace.decToBin(input1); 

cout << "\n"; 

necklace.printSequence(); 

necklace.findDivisors(length); 


cout << "\n\n"; 
cout << "Period = " << necklace.findSequencePeriod() << endl; 


cout << "\n\nDo you want to find the period of another sequence? "; 
cin >> answer; 


} 


cout << "\n\nDo you want to find the necklace of a sequence? "; 
cin >> answer; 
while (answer == 'Y' || answer == 'y’) 


{ 


cout << "\nPlease enter a decimal to generate a binary sequence: "; 
cin >> input1; 

necklace.decToBin(input1 ); 

cout << "\nSequence prior to shifting: "; 
necklace.printSequence(); 


cout << "\n"; 
necklace.shiftSequenceToNecklace(); 


cout << "\nNecklace: " << endl; 
cout << "\n"; 


necklace.printSequence(); 
cout << "\n\nDo you want to find the necklace of another sequence? "; 


cin >> answer; 


} 


cout << "\n\nDo you want to see the distances between the shifting ones signposts? "; 
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cin >> answer; 

if (answer == 'Y' Il answer == 'y') 

{ 
cout << "\n"; 
necklace.displayShiftingOnesDistances(); 


} 


cout << "\n\nDo you want to see the deBruijn section surrounding your necklace input? 


cin >> answer; 
while (answer == 'Y' Il answer == 'y’) 
{ 
cout << "\nPlease enter a decimal between 0 and " << 
necklace.generatePower2(length)-1 <<" to generate a binary sequence: "; 
cin >> input3; 
necklace.decToBin(input3); 
cout << "\nSequence prior to shifting: "; 
necklace.printSequence(); 
cout << "\n"; 
necklace.shiftSequenceToNecklace(); 
cout << "\nNecklace: " << endl; 
cout << "\n"; 
necklace.printSequence(); 
input3 = necklace.binToDec(); 
cout << "\n\ndecimal associated with necklace = " << input3 << endl; 
cout << "\n\nPlease enter another decimal greater than " << input3 << ":"; 
cin >> input1; 
necklace.decToBin(input1 ); 
cout << "\nSequence prior to shifting: "; 
necklace.printSequence(); 
cout << "\n"; 
necklace.shiftSequenceToNecklace(); 
cout << "\nNecklace: " << endl; 
cout << "\n"; 
necklace.printSequence(); 
inputl = necklace.binToDec(); 
cout << "\n\ndecimal associated with necklace = " << inputl << endl; 
cout << "\n\nPlease enter another decimal less than " << input3 << ":"; 
cin >> input2; 
necklace.decToBin(input2); 
cout << "\nSequence prior to shifting: "; 
necklace.printSequence(); 
cout << "\n"; 
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necklace.shiftSequenceToNecklace(); 
cout << "\nNecklace: " << endl; 

cout << "\n"; 
necklace.printSequence(); 


cout << "\n\nRunning Necklace Algorithm..." << endl; 
necklace.runNecklaceAlgorithm(input1, input2, input3, false); 


cout << "\nFirst necklace: "<< endl; 
necklace.displayFirstNecklace(); 


cout << "\nMiddle necklace: " << endl; 
necklace.displayMiddleNecklace(); 


cout << "\nLast necklace: " << endl; 
necklace.displayLastNecklace(); 


necklace.createDeBruijnSectionQ); 


cout << "\nDeBrujjn Section..." << endl; 
cout << "\n"; 
necklace.displaydeBruijnSectionQ); 


cout << "\n\nDo you want to search the deBruijn section for a particular " << length 
<< "-tuple? "; 
cin >> answer; 
while (answer == 'Y' Il answer == 'y’) 
{ 
int pos = 0; 
int per = 0; 
cout << "\nPlease enter a decimal value between " << inputl <<" and " << input2 
ee Ne 
cin >> input3; 
necklace.getInputDec(input3); 
cout << "\nSequence prior to shifting: "; 
necklace.printSequence(); 
cout << "\n"; 
necklace.shiftSequenceToNecklace(); 
cout << "\nNecklace: " << endl; 
cout << "\n"; 
necklace.printSequence(); 
input3 = necklace.binToDec(); 
cout << "\n\ndecimal associated with necklace = " << input3 << endl; 


necklace.runNecklaceAlgorithm(input1, input2, input3, true); 
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cout << "\nFirst necklace in deBruijn section: " << endl; 
necklace.displayFirstNecklace(); 
necklace.createInputString(input3); 
cout << "\nBinary string associated with decimal input: "; 
necklace.printSequence(); 
cout << "\n"; 
necklace.createDeBruijnSection(); 
pos = necklace.searchDeBruijnString(); 
cout << "\n\nPositon associated with binary string input = " << pos << endl; 
cout <<"\nDo you want to search for another sequence? " << end]; 
cin >> answer; 
} 
cout << "\n\nDo you want to see another deBruijn section surrounding your necklace 
input? "; 
cin >> answer; 


} 


cout << "\n\nGoodbye.." << endl; 
IOThesis::pauseBeforeExit(EXIT_NORMALLY); 


return 0; 
} // end Thesismain.cpp 


3. IO Thesis CPP File 
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// 1OThesis.cpp 

// 

// LT John Ortiz 

// 

// Thesis: Generating Necklaces 

// Operating Environment: Windows XP Home 
// Compiler: Visual Studio .NET 

// Date: 08 December 2005 


// Description: This program verifies that user inputs are valid 
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#include <iostream> 

#include <cstdlib> 

#include <iomanip> 

#include <cctype> 

#include <string> 

#include <sstream> 

#include "IOThesis.h" 

#include "FileOpeningException.h" 


using namespace std; 
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using std::exit; 
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// Method displayMessage 

// Return value none 

// Parameters none 


// Purpose prompts user for dimension input 
J [78 A AR A A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 2 2 ee he he fe fe fe he he ie Ie ie ie ie 2S ik 


void IOThesis::displayMessage(string message) 
{ 


cout << message << end]; 


} //end displayMessage 
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// Method openOutputFile() 

// Return value 

// Parameters: const char *FILE_NAME the name of the file to open for writing 


// Purpose: dynamically create a new ofstream and returns a pointer if successful 
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ofstream* IOThesis::openOutputFile(const char *FILE_NAME) 


{ 
ofstream *outputFileptr = NULL; 
try 
{ 
outputFileptr = new ofstream(FILE_NAME,ios::out); 
if (!(*outputFileptr)) // unable to open object outputFile of type ofstream 


{ 
throw FileOpeningException(); 


} 


else 


{ 


return outputFileptr; //returns pointer to ofstream if file was opened successfully 
} 


} 
catch(FileOpeningException &) 
{ 

throw; 


} 
}// end openOutputFile 
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// Method openInputFile 

// Return value ifstream* 

// Parameters const char *FILE_NAME, the file to open for reading 
// Purpose 
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ifstream* IOThesis::openInputFile(const char *FILE_NAME) 
{ 
ifstream *inputFileptr = NULL; 
try 
{ 
inputFileptr = new ifstream (FILE_NAME,ios::in); 
if (!(*inputFileptr)) // unable to open object outputFile of type ofstream 
{ 
throw FileOpeningException(); 


} 


else 


{ 
return inputFileptr; //returns pointer to ofstream if file was opened successfully 
} 
} 
catch(FileOpeningException &) 
{ 
throw; 


} 
} // end openInputFile 
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// Method inputInt 

// Return value int 

// Parameters none, function will ask for input via I/O 


// Purpose take string input and determine if valid for an int (integer can be any value 
less than MAXINT) 
// use inputFilter = [OThesis::inputInt(); vice cin >> inputFilter 
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int IOThesis::inputInt() 


{ 
const int SIZE = 12; 
char ch[SIZE]; // temporary buffer for input characters 
char ch2[SIZE]; // buffer for sstream 
string token; 
char *tokenPtr; 
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int size = 0; 

int integer = 0; 
int count = 0; 
bool state = false; 
stringstream ss; 


while(state == false) // This loop prompts a user to enter in only one integer that is less 
than 12 characters long 
{ 
int q = 0; 
int total=0; 
char element = cin.get(); 
while(element != ‘\n')// reads elements into a character array 
{ 
ch[q]= element; 
total++; 
if (total == SIZE) // checks if total exceeds 12 characters 
{ 
cout << "You have reached the maximum number of characters allowed" << 
endl; 
cout << "Please enter another integer."; 
state = false; 
cin.clear(); 
break; 
} 


if (cin.eof()) // checks if ctrl-z in input as a character 


{ 


cout << "\nInvalid character."; 
cin.clear(); //resets input stream so it can keep processing 


cin.putback(‘\n'); 
state = !IOThesis::clearInputBufferQ); 
break; 
} 


if (q == 0 && (ch[q] == '+'ll ch[q] =='-' )) // checks if first character is plus or 
minus 
{ 
element = cin.get(); // gets the next element 
q++; 
continue; 


} 


else if(!isdigit(ch[q])) //checks if I have something other than a digit 
{ 
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cout << ch[q] << " is not an integer."<< endl; 
state = IOThesis::clearInputBuffer(); 
state = false; 
cin.clear(); 
break; 
} 
else 
{ 
state = true; 
element = cin.get(); // gets the next element 
q++; 
} 
} // end while 


tokenPtr = strtok (ch," "); //begin tokenization of string 
int numtok = 0; // need to initialize through every pass of outer while loop 
while (tokenPtr != NULL) // loop exits if there is one or more token detected 
{ 
tokenPtr = strtok (NULL, " "); //get next token 
numtok ++; 
} 7 — end while 


if (numtok != 1) // tests if there is more than or no inputs in the string 


{ 


cout << endl; 
cout << "You have too many inputs." << endl; 
cin.putback(‘\n'); 
state = IOThesis::clearInputBuffer(); 
state = false; 
cin.putback(‘\n’); 
} // end if 
else 
{ 
state = true; 
count = q; 
} 


} //ends outer while loop 
if (state == true) 
for (int m=0;m < count; m++) 
ch2[m] = ch[m]; // only assigns current "good" input to sstream 
} 


ss << ch2; //takes input string and sends it to sstream 
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Ss >> integer; 
} // end if 
return integer; 


} // end inputInt 


J [787 AR A A A 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 2 ee fe he fe he fe he fe ik 2c ie ie 2s Ie iS 


// Method inputIntWithLimits 

// Return value int 

//Parameters (string message, int lower, int upper) 

// Purpose calls input int for a valid int, then checks limits 


// used similiarly to inputInt 
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int IOThesis::inputIntWithLimits(string message, int lower, int upper) 
{ 
cout << message << endl; 
int condition = IOThesis::inputInt(); 
while ( condition < lower II condition > upper) 
{ 
cout << message << endl; 
condition = IOThesis::inputIntQ); 
}// end while 


return condition; 


}// end inputIntWithLimits 
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// Method clearInputBuffer 
// Return value bool TRUE if cin buffer has no additional stuff in it 
//Parameters none 


// Purpose utility function for processing inputs 
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bool IOThesis::clearInputBuffer() 
{ 
bool boolstate = false; 
char somechar = cin.get(); 
while(somechar!= '‘\n') //loop enables the input buffer to be flushed 
{ 
if (isspace(somechar)) 


{ 
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boolstate = true; 


} 
else 
{ 
boolstate = false; 
} 


somechar = cin.get(); // keeps writing over somechar's memory space with current 
character 
} // end while 


return boolstate; 


} //end clearInputBuffer 
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// Method pauseBeforeExit 

// Return value none 

// Parameters int condition (EXIT_WITH_ERROR = 1, global defined in IOThesis.h) 
// (EXIT_NORMALLY = 0, global defined in IOThesis.h) 


// Purpose keeps display window open long enough to see results 
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void IOThesis::pauseBeforeExit(const int condition) 
{ 
cout << "\n\n" << endl; 
char ch; 
ch = cin.get(); //need to flush the input buffer or it may still hold the last 
// carriage return and blow us out of the program early 
cout << "\n\tPress <enter> or <return> to exit...\n"; 
ch = cin.get(); 
exit(condition); 


} //end pauseBeforeExit 


// end of file IOThesis.cpp 
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